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r-| ' ABSTRACT 

Or 

Aims. We present results from a deep spectral analysis of all the Swift observations of Mrk 421 from April 2006 to July 2006, when it 
reached its largest X-ray flux recorded until 2006. The peak flux was about 85 milli-Crab in the 2.0-10.0 keV band, with the peak energy (E p ) 
of the spectral energy distribution (SED) laying often at energies larger than 10 keV. We study the trends among spectral parameters, and their 
r ^2 ' physical insights in order to understand the underlying acceleration and emission mechanisms. 

1 Methods. We performed spectral analysis of the Swift observations investigating the trends of the spectral parameters in terms of acceleration 

and energetic features phenomenologically linked to the SSC model parameters, predicting their effects in the y-ray band, in particular the 
spectral shape expected in the Fermi Gamma-ray Space Telescope-LAr band. 

Results. We confirm that the X-ray spectrum is well described by a log-parabolic distribution close to E p , with the peak flux of the SED (S p ) 
being correlated with E p , and E p anti-correlated with the curvature parameter (b). The spectral evolution in the Hardness-ratio-flux plane shows 
both clock-wise and counter-clock-wise patterns. During the most energetic flares the UV-to-soft-X-ray spectral shape requires an electron 
distribution spectral index s 2.3. 

Conclusions. Present analysis shows that the UV-to-X-ray emission from Mrk 421 is likely to be originated by a population of electrons 
that is actually curved, with a low energy power-law tail. The observed spectral curvature is consistent both with stochastic acceleration or 
energy dependent acceleration probability mechanisms, whereas the power-law slope form XRT-UVOT data is very close to that inferred from 
the GRBs X-ray afterglow and in agreement with the universal first-order relativistic shock acceleration models. This scenario hints that the 
magnetic turbulence may play a twofold role: spatial diffusion relevant to the first order process and momentum diffusion relevant to the second 
order process. 
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1. Introduction 

BL Lac objects are Active Galactic Nuclei (AGNs) charac- 
terized by a polarised and highly variable nonthermal contin- 
uum emission extending from radio to y-rays. In the most ac- 
cepted scenario this radiation is produced within a relativis- 
tic jet originated by the central engine and pointing close to 
our line of sight. The relativistic outflow has a typical bulk 
Lorentz factor r w 10, hence the emitted fluxes, observed at 
an angle 6, are subjected to the effects of a beaming factor 
8 = 1/(T(1 -f3 cos 6)). 



The Spectral Energy Distribution (SED) of these objects 
has a typical two-bump shape. According to current models, 
the lower-frequency bump is interpreted as synchrotron emis- 
sion from highly relativistic electrons with Lorentz factors y 
in excess of 10 2 . This component peaks at frequencies ranging 
from the IR to the X-ra y band. The actual position of this peak 
has been suggested by Padovani & Giommi (1 1995b as a marker 
for a classification; they define LBL (Low energy peaked 
BL Lac) the objects with the first bump in the IR-to-optical 
band, and HBL (High energy peaked BL Lac) those peak- 
ing in the UV-X-ray band. According to the Synchrotron Self 
Compton (SSC) emission mechanism, the higher-frequency 
bump is to be attributed to inverse Compton scattering of syn- 
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Table 1. Swift observation journal and exposures of Mrk 421. 
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(*) Pointings triggered by the BAT instrument as GRBs. 



chrotron photons by the same population o f relativistic elec- 
trons that produce the synch rotron emission ( Jones et aDll974 ; 
Ghisellini & Maraschill 1989b . 

With its redshift z = 0.031, Mrk 421 is among the closest 
and best studied HBL. In fact, it is one of brightest BL Lac 
objects in th e UV and in th e X-ray bands, observed in y rays 
by EGRET (ILin et al.lll992l) : it was also the first extragalactic 
source detected at Te V energies in the range 0.5- 1 .5 Te V by the 
Whipple telescopes dPunch et alJll992l;|Petrv et alJll996l) . 

The source is classified as HBL because its synchrotron 
emission peak ranges from a fraction of a keV to several 
keV. Its flux changes go along with strong spectral varia- 
tions (Fos sati et al.l2 000a; M assaro et al . 2004) and the spectral 
shape generally exhibits a marked curvature, well described by 



log-parabolic distributio n with a mildly cu rved and symmetric 



2007 



;-pa 
b). 



a log-p arabolic model dMassaro et al.l l2004t iTramacere et al 



In spring/summer 2006 Mrk 421 reached its largest X-ray 
flux recorded until that time. The peak flux was about 85 
milli-Crab in the 2.0-10.0 keV band, with the peak energy of 
the spectral energy distribution (SED) often lying at energies 
larger than 10 keV. 

In this paper (Paper I) we present data collected from Swift ob- 
servations performed during this very intense flaring period. 
We aim to study the evolution of the spectral parameters as a 
function of the flaring activity, and the correlations among the 
spectral parameters. This gives a phenomenological picture 
of the physical mechanism driving the observed patterns. In 
Paper II we will frame th is scenario in the theoretical context 
of stochastic acceleration ( Tramacerel2009l) . 



In the phenomenological context of jets in HBLs the spec- 
tral curvature is relevant for understanding of both radiative 
and acceleration mechanisms. 

Many works in the literature have shown that the X-ray spec- 
tral shape of Mrk 421 is actually curved and described by a 



spectral shape dFossati et al.l l2000rx iTanihata et al 
Massa roet al.l 120041 : ITramacere et al.l l2007ah . iMassaro et al 



2004 



( 2004) gave an interpretation of this feature in the framework 
of energy dependent acceleration efficiency that naturally leads 
to log-parabolic spectral distribu tions with a possible power- 
law tail at lower energies. iKardashevI {[962) showed that a 
log-parabolic distribution results from a stochastic acceleration 
scenario with a mono-energetic or qua si-m ono energetic 



Katarzvhski et al. (2006) and Giebels et al 



particle injection 
(2007) used a relativistic Maxwellian electron distribution 



resulting from a stochastic acceleration process, to describe the 
X-ray /Te V emission of Mrk 501 and Mrk 421 respectively. 
Recently. IStawarz & Petrosian (2008) showed that a distribu- 
tion with similar spectral properties can be obtained as a steady 
state energy spectra of particles undergoing momentum dif- 
fusion due to resonant interactions with turbulent MHD modes. 



Tramacere et al. (2007b) suggested that the connection of 



the X-ray curvature with that in the emitting particles and its 
evolution with the source state, could be investigated as a test 
for the prediction of the scenarios li sted above. In particular, 
both stochastic dKardashevI 1 19621) and energy-dependent 



acceleration mechanisms predict an anticorrelation between 
the curvature and the SED peak energy (E p ). 
Moreover, the pattern shown by the peak height (S p ) of the 
SED as E p moves, can trace the evolution of the parameters 
characterizing the energetics of the synchrotron emission, in 
particular the average particle energy and the number density 
of the emitting particles. 



A crucial issue is understanding whether during the most 
violent flares the shape of the X-ray spectrum also can be 
described by a single log-parabolic spectral distribution. In 
fact, a typical X-ray detector shows only a slice (usually up 
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Fig. 1. From top to bottom: a) light curve of the flux in the band 0.2 - 10.0 keV. b) light curves for three different bands soft 
(0.2 - 3.0 keV, red points) medium (3.0 - 5.0 keV, green points) and hard (5.0 - 10.0 keV, blue points), normalized to their 
maximum value . c) the evolution of the hardness ratio (HR) evaluated as the ratio of the 2.0 - 10.0 keV band to the 0.2 - 2.0 keV 
band, d) light curve from the Swift-UVOT instrument, different colors refers to different filters (V=brown, U=black, UVM2=red, 
UVWl=blue, UVW2=yellow). [See the electronic edition of the Journal for a color version of this figure] 



to two decades) of the overall emission from the observed 
object. During strong flares involving large variation of the 
SED peak energy, it is possible to understand if the electron 
distribution is curved even far from the peak energy. Moreover, 
we can capitalise on the unique opportunity given by Swift to 
perform simultaneous UV-to-X-ray observations, extending 
the spectral window from about 10 15 Hz up to 10 19 Hz. The 
presence of a power-law tail at low photon energies and its 
slope can provide information about the low-energy tail of the 
underlying electron distribution as well as on the acceleration 
mechanism generating such a spectral shape. 



As final task, building on the phenomenological re- 
sults from the present analysis, we can model the SED 
of Mrk 421 within the synchrotron-self-Compton (SSC) 
scenario, predicting the possible spectral behaviour at y-ray 
energies. In particular we can relate the typical spectral shape 
of the UV-to-soft-X-ray emission to that expected in the 
MeV/GeV band covered by the Fermi Gamma-ray Space 
Telescope-LAr instrument. 

This paper is organized as follows. In Sect. 2 we present 
our data set and the procedure used to reduce the Swift data. In 
Sect. 3 we report results from the analysis of the fast variability 
of the source. In section 4 we present spectral analysis results. 
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In section 5 we analyse the spectral evolution of the source. 
In Sect. 6 we study the trends among the spectral parameters 
and compare these results with previous studies directing our 
attention to the link between such trends and expectations from 
different scenarios for acceleration mechanism. In Sect. 7 we 
study the connection between the UVOT and the XRT spec- 
tra, showing the relevance of the derived spectral shape in the 
context of the first order acceleration processes. In sect. 8 we 
model the SED within the SSC framework, focusing on the 
phenomenological interpretation of the data. In Sect. 9 we dis- 
cuss the results, and in Sect. 10 we draw our overall conclu- 
sions. 

2. Swift Observations and data analysis 

We present results from temporally-resolved spectral analyses 
of 15 Swift simultaneous observations in the UV/X-ray band 
performed between April and July 2006, when the source was 
so bright that the high-energy instrument BAT automatically 
triggered pointings at it three times assuming as if it were 
a Gamma Ray Burst. The log of UV/X-ray observations is 
reported in Tab. [T] 

In Fig. [T] we report the XRT light curves obtained from the 
single-orbit spectra. In the first panel from the top we show 
the light curve of the flux obtained integrating the model from 
Eq. |2] (see section 4.2) between 0.2 and 10.0 keV, according 
to the parameters and parameters errors reported in Tab [2] 
Fluxes in Tab. 2 refer to the 0.3-10.0 keV interval because 
the XRT response function is calibrated only in that range. 
We extrapolated the flux to the 0.2-10.0 keV band to make 
easier comparisons with data from other X-ray telescopes that 
usually are given in the 0.2-10.0 keV band. 
The second panel from the top shows the light curve for three 
different bands: soft (0.2 - 3.0 keV) medium (3.0 - 5.0 keV) 
and hard (5.0 - 10.0 keV), normalized to their maximum value. 
In the third panel from the top, we report the evolution of the 
hardness ratio (HR) evaluated as the ratio of the 2.0 - 10.0 keV 
band to the 0.2 - 2.0 keV band. 

The bottom panel of Fig. Q] shows the light curve obtained 
from the Swift-UVOT observations. 



of the source (40-80 cts/s) . 

We selected photons with grades in the range 0-2 for WT 
mode; we also used default screening parameters to produce 
level 2 cleaned event files. To take advantage of the good statis- 
tics offered by such large numbers of events we decided to 
make a high temporal resolved analysis, extracting spectra for 
each Swift orbit. We rejected only two out of the 174 obtained 
spectra, because of a strongly biased exposure. The resulting 
XRT database presented in this work therefore includes 172 
time intervals. 



2.2. Swift-UVOT data analysis 



The Swift UV and Optical Telescope (UVOT lR^min3l2005b 
observations included in this paper have exposures in all fil- 
ters except for the White one. Photometry of the source was 
performed using the standard UVOT software developed and 
distributed within the HEAsoft 6.4 package. Counts were ex- 
tracted from an aperture of 5" radius for all single exposures 
within an observation and for all filters, while the background 
was carefully estimated in few ways. In alm ost all observations 
the source is on the "ghost wings" ( O 20061) from the nearby 
star 5 1 UMa, so we estimated the background from a circular 
aperture of 15" radius off the source but on the wings, exclud- 
ing stray light and support shadows. These background values 
were compared to those obtained using a region outside the 
wings, resulting in differences in most cases within the errors. 

We checked astrometry of each exposure, verifying the 
aperture positioning. Count rates were then converted to fluxes 
using the standard zero points. We discarded some exposure for 
which the count rate was near the limit of acceptability for the 
"coincidence loss" correction factor included in the CALDB 
(-90 cts s- 1 ). 

The fluxes were then de-reddened using a value for E(B-V) 



of 0.014 mag (Schl egel etal.ll 1998b with A A /E(B - V) ra 



tios calculated for UVOT filters (for the latest effective wave- 
lengt hs) using the mea n Galactic interstellar extinction curve 
from iFitzpatrickl dl999h . 



2.1. Swift-XRT data analysis 



2.3. Swift-BAT data analysis 



All the data were reduced using the XRT DAS software 
(version v2.2.0) developed at the ASI Science Data Center 
(ASDC) and distributed within the HEAsoft 6.4 package by 
the NASA High Energy Astrophysics Archive Research Center 
(HEASARC). 

The operational mode of XRT is automatically controlled 
by the on-board software that uses the appropriate CCD read- 
out mode to reduce or eliminate the effects of photon pile-up. 
When the target count-rate is higher than wl cts/s the system is 
normally operated in Windowed Timing (WT) mode whereas, 
the Photon Counting (PC) mode is used for fainter sources 



(see iBurrows et alJl2005t iHill et al.ll2004 for more details on 



XRT observing modes). The observations presented here were 
all performed in WT mode, owing to the extremely high state 



We analysed the data that automatically triggered the BAT in- 
strument (labeled with *) in Tab. [TJ. To reduce the data we 
followed the instructions reported in the BAT analysis threads 
http : //swift . gsfc . nasa . gov/docs/swi ft/analysis/ 
threads/batspectrumthread.html. We first gener- 
ated a detector plane image (DPI) and checked for noisy 
pixels using the bathotpix task. In order to properly 
subtract the background we generated a weighting mask 
using the batmaskwtevt task. The spectrum was ex- 
tracted using the batbinevt task. We applied geometric 
corrections (batupdatephakw), and systematic error correc- 
tion(batphasyserr). As final task we generated a response 
matrix (batdrmgen). The spectrum was rebinned by providing 
an external bin edges file to batbinevt. 
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3. Fast source variability. 

Temporal variability is one of the most interesting features 
characterizing HBLs. In the X-ray band typical time scales for 
flux changes decrease to the order of hour or minutes. The iden- 
tification of the order of the fastest significant flux change al- 
lows to estimate the source size, given by the well-known rela- 
tion: 

c At 6 



R < 



(1) 



1 +z ' 

where c is the speed of the light, 6 the beaming factor, and z 
the cosmological redshift. We investigated the intra-orbit light 
curves from our data set and found a significant flux variation 
in the first orbit on 04/26/2006 (ObsID 00030352006) pointing 
(see Fig.|2]i. The light curve extracted from the 0.3 - 10.0 keV 
band shows an increasing trend spanning about 1200 seconds 
with a variation of about 35%. 



Bin time: 50.00 s 



analysing TeV data from the CAT telescope (R < 5.4 x 10 13 r5 
cm). 



4. Spectral analysis 

4.1. XRT spectral analysis 

We find, for most of the spectra in our sample, systematic 
deviations (see Fig. |3]l in the residuals obtained when fitting 
the data by means of a power-law spectrum with Ah fixed 
at the Galactic value. Such a behaviour heuristically suggests 
that the spectra are intr insically curved. Th i s was already 



known in the literature dFossati et"aL 2000bl iTanihata et aL 
120041 iMassaro et al . 120041; iTramacere et al.l2007allbl) . All these 
authors agreed that when the spectral shape of Mrk 421 is 
curved, describing the curvature only in terms of absorption 
not only would require a column density m uch higher than the 
Galac tic value Nh = 1.61 x 10 20 ctrT 2 dLockman & Savage 
1995b . but also would yield in any case unacceptable fits with 
very high x 2 r - Moreover, high resolution images of the host 
early-type galaxy of Mrk 421 do not show in the brightness 
profile any evide nce of large amounts of absorbing material 
(lUrrv et al. I l2000l) . Based on these phenomenological results, 
we performed the spectral analysis fixing the Ah absorbing 
column densities to the Galactic values and using the following 
log-parabolic spectral law (LP): 



F(E) = KE 



-(a+b log(£)) 



ph cm 



- 2 s- 1 keV- 1 , 



(2) 



where a is the photon index at 1 keV and b measures the 
spectral curvature. 



2000 



2200 



2400 



2600 
Time (s) 



2800 



3000 



3200 



Start Time 53851 (MJD) 3:32:32:174 Stop Time 53851 (MJD) 3:53:22:174 

Fig. 2. Light curve binned with 50 s intervals, from the first or- 
bit of the 26/04/2006 pointing. Significant variability (^ 35%) 
is detected over intra-orbit time scale (1200 s). 

According to Eq. Q] a 1200 s time scale implies 
R < 3.6 x 10 13 r5 cm. Assuming a beaming factor of the 
order of 10 we get R 4 x 10 14 cm, which indicates a quite 
compact em it ting re gion. 



Lichti et al.l d2008l) analysed the temporal variability of 
Mrk 421 using observations performed by the INTEGRAL- 
ISGRI instrument in the 40 - 100 keV band, and overlapping 
our data set during the June pointings. The fastest variability 
observed in the INTEGRAL-ISGRI light curves, estimated 
by fitting the data with a rise-time law of the form a ■ e'^'", 
gives R < 3 x 10 14 r5 cm. We perform the same analysis for 
the 04/26/2006 light curve and we find t 4.100 s which 
implies R < 1 x 10 14 <5 cm. We use this time scale to constrain 
the emitting region size in the following analyses (Sect. 8). 
Time scales si milar to that in the 04 /26/2006 pointing were also 
~~ d200 lb in ASCA data , about 5 ks 



observed by Tanih ata et al. 



The SED peak energy (E p ) and the SED height (S p ) can 
be derived easily from Eq. [2] but in this case they suffer in in- 
trinsic analytical correlation. This bias can be removed using 
an equivalent functional relationship that is a log-parabola ex- 
pressed in terms of E p , S p and b (LPEP): 



erg cm 2 s 1 , 



(3) 



(«<1.5x 10 l4 5 cm), and by Gie bels et all (120071) (about 2 ks) 



S(E) = (1.60 x 10~ 9 ) S p 1(T* Cos^)) 2 

where S p = E p F(E p ) and E p are estimated during the fit, and 
the numerical constant is simply the energy conversion factor 
from keV to erg. 



4.2. Orbit resolved analysis 

Because of the very bright state of the source, we were able 
to extract spectra for each orbit, for a total of 172 spectra. 
A motivation to perform an orbit resolved analysis is the 
strong variability of the source during these pointings. In fact, 
integrating spectra over time scales much longer than the 
typical variability leads to misleading results in the estimates 
of the curvature, E„, and S p . 

The results of the spectral analysis are reported in Tab. 
[2] (rejected spectra are labeled with (*)), where all statisti- 
cal errors refer to the 68% confidence level (one Gaussian 
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0.5 1 2 5 0.5 1 2 5 

channel energy (keV) channel energy (keV) 



Fig. 3. Spectrum from the first orbit of the ObsID 00030352013 performed on 06/22/2006. Left Panel: the systematic deviations 
on both sides of the residuals from a best fit with a power-law with Galactic Nh show the need of intrinsic curvature. Right 
Panel: the deviations disappear with the log-parabolic model with Galactic Nh ■ The xl decreases from 1.60 with 246 d.o.f. 
(power-law) to 1.19 with 245 d.o.f. (log-parabola); the F-test statistics clearly favours the curved model. 



standard deviation). The second thorough forth columns in 
Tab. [2] report the best fit estimates for the model in Eq. [2] 
The fifth column reports the value of the SED peak analyt- 
ically estimated from Eq. [2] according to the best fit results 
(E p * ). The sixth and seventh columns report the E p and 
S p best fit estimates using as fitting model Eq. [3] In the 
eighth column we report the flux in the 0.3-10.0 keV band, 
evaluated by X-spec integrating the Eq. [2] model. In the last 
column we report the reduced^ 2 statistics for the fit with Eq.|2] 

The SED peak energy was often difficult to estimate. The 
reason was that during this particular high brightness state, the 
spectra were in some cases very hard, with a photon index 
a =i [1.6-1.7] and and having a low spectral curvature, im- 
plying a peak energy placed far out of the XRT energy band. 
In order to test the robustness of the E p estimate, we first de- 
rived the peak energy from the spectral parameters of Eq. [2] 
{E p *). Then we fitted the spectra using Eq.[3] setting the initial 
curvature value to that returned from the fit with Eq. [2] In or- 
der to test the stability of the results we adopted the following 
criteria: 



1. The value of E p statistically significant. Given the asym- 
metric uncertainties we define a- Ep the half 2 sigma confi- 
dence level, and require E p /(Te p < 1. 

2. E p * consistent with E p at within one sigma. 

We report in Tab.|2]the estimates of E p satisfying this criterion, 
in the other cases we report only the lower limit of E p *. The 
estimates of E p * > 100 keV obviously are not statistically 
robust, meaning that the actual energy peak may be in excess 
of 100 keV, but we could not give a robust estimate. 
All the spectra for which the stability conditions were satisfied 
returned values of E p < 20 keV. 




Fig. 4. Joint XRT BAT spectral analysis from the 04/22/2006 
pointing. 



4.3. Orbit merged analysis 

The orbit resolved spectral analysis has the great advantage to 
follow accurately the strong variability of the source, but the 
E p estimates suffer from large uncertainties. In any case, based 
on the spectral/flux pattern traced by the previous analysis, we 
can identify all the orbits indicating essentially the same spec- 
tral/flux states. We can use these intervals to perform an orbit- 
merged spectral analysis, to get smaller uncertainties on the E p 
value, and without integrating the source over periods exhibit- 
ing large changes. 

The result of this analysis is reported in Tab. [3] In this analysis 
when E p and E p * can be determined the typical uncertainties 
are smaller. We use the values from this table to perform the 
E p - S p and E p - b trend analysis in the following, selecting 
only the data with significant E p estimates. Also in this case, 
all the spectra with E p well constrained have E p < 20 keV. 
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2.0x10 2.5x10 3.0x10 3.5x10 

flux (erg cm s ) 

Fig. 5. Scatter plot of the HR vs. the flux in the 0.2-10keV band 
(HR is evaluated as the ratio of the 2.0 - 10.0 keV band to the 
0.2 - 2.0 keV band). The points in the shaded area are almost 
all from the period 06/23/2006 to 06/27/2006 and 07/15/2006. 

4.4. Joint XRT-BAT spectral analysis 

For the three observations that have simultaneous Swift- 
XRT and Swift-BAT data (with the BAT in automatic trigger 
mode, see Tab. [TJ we performed a joint XRT-BAT spectral 
analysis. The results are reported in Table [3] (lines labeled 
X+B). The spectral curvature resulting from the joint analysis 
is slightly larger with respect to to the case of only XRT data, 
but is always consistent within one sigma. As already dis- 
cussed, the estimation of E p is more subtle. We will use 
the superscript X+B to refer to joint XRT and BAT analysis 
hereafter. 

For the 07/15/2006 pointing, the two curvatures are 
b = 0.17 + 0.02 and b x+B = 0.20 + 0.02, and the two 
values of E p consistent within one sigma are E p = 11 +~ ke V 
and E x+B = 8^ keV. For the case of the 04/22/2006 pointing 

20+™ keV 
0.11 + 0.02 and 



(Fig. |4]i, the peak energies resulting are E XJ[ 



and E„ 



26+1 9 and the curvatures b 



b x+B — 0.12 + 0.02. 
For the 06/23/2006 pointing, the XRT curvature is estimated 
0.08 + 0.03 to be compared to the value of b x+B = 0.13 + 0.02. 
In this case the smaller value of the curvature, as discussed in 
the previous section, makes more difficult the estimate of the 
E p value. In fact, XRT data are not able to locate the actual 
value of the peak energy, and the estimate from joint analysis 



is E x+B 



34+22 ke y 



The results from joint XRT-BAT analysis confirm that E p 
values estimated from XRT data, at energies larger than about 
20 keV, are potentially biased. The selection applied in the 
present analysis should warrantee that the bias on the trends 
among E p , b and S p will be as small as possible. 

5. Spectral evolution 

In this section we follow the classical approach of the hardness 
ratio (HR) analysis, to investigate cooling and/or acceleration 
features. The HR is evaluated as the ratio of the 2.0 - 10.0 keV 
band to the 0.2 - 2.0 keV band. 



Analysis of the evolution of the HR shows that the correlations 
between variations in the soft and hard bands results in 
modulation of the HRs, with the spectra harder when the 
source is brighter and softer when weaker. This trend is clearly 
visible in Fig. [5] A significant scatter in the points hints that 
the dynamic among different flares is quite different due to 
different underlying physical conditions. 
Moreover, we note that all but 3 of the points in the shaded 
area belong to the period from 06/23/2006 to 06/27/2006 and 
to the 07/15/2006 pointing, when the source was brighter than 
in the previous pointings. 

A deeper understanding of the spectral dynamics can be 
achieved from looking at the hysteresis patterns of the single 
flares in the a-flux or HR-fixix plane, where a is the photon in- 
dex at 1 keV (see Tab. [2j. 

The time scales relevant to understanding the patterns in the 
a-flux plane are: the injection (r,„y), the escape (j esc ), the cool- 
ing (t coo i), the acc eleration (T n (]r ) and the light crossing (r eross ) 
time. According to Kir k et al. ( 1998b . the loops are expected to 
be clockwise (CW) and with soft lag when the flare is observed 
at frequencies where the higher energy variability is faster than 
at lower energy, (as in the case of synchrotron cooling). On 
the contrary, when observed at frequencies for which the accel- 
eration and cooling time scale are almost equal, the loops are 
expected to be counter-clock-wise (CCW) with a possible hard 
lag. 

We investigated carefully all the patterns from our data-set, and 
in Fig.|6]and Fig.|7]we report the 5 flares showing a clear CW or 
CCW loop. The left panels show the evolution of a as a func- 
tion of the flux, and the right panels show the light curves in 
three different bands. The fluxes are normalized to the maxi- 
mum value reached in that particular time interval. As a general 
result, we did not find any significant plateau among the light 
curves of these flares, meaning that t,„ ; - 
the behaviour of the flares individually: 



< T, 



We describe 



Flare 1) The flare has a CCW pattern; it begins with a soften- 
ing of the spectrum with a flux almost steady (A) probably 
reflecting the cooling from the previous flaring episode, fol- 
lowed by a flux increase with a spectral hardening (B). 

Flare 2) The flare has two patterns: one CW and one CCW. It 
starts with a decreasing flux and a spectral softening (A) ( 
Tcooi ^ Tesc)- Figure [6] (left panel) shows that the flux in the 
soft band is still increasing when the fluxes in the medium 
and hard bands are already decreasing. When the spectrum 
starts to get harder (B) the flux is still decreasing, this may 
hint that we are seeing the propagation of the new injec- 
tion starting from the hard band, with the soft band still 
decreasing. Then, the flux starts increasing with the spec- 
trum hardening (C). When the flare starts to decay (D) the 
pattern switches to CCW and we see a decrease of the flux 
with a spectral softening. 

Flare 3) The flare has a CCW pattern, starts with a flux de- 
crease and a spectral softening (A), then the flux increases 
and the spectrum hardens (B). In the last part (C), the flux 
decreases and the spectrum still hardens. As in the previous 
case this may hint the start of a new hard flaring component. 
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Fig. 6. Left panels: Spectral patterns in the a-flux plane showing clockwise and counterclockwise trends. Right panels: Swift- 
XRT light curves at three different bands: soft (0.2 - 3.0 keV, red), medium (3.0 - 5.0 keV, green), and hard (5.0 - 10.0 keV, red) 
. Capital letters guide the reader during the comment in the text. In both left and right panels, the fluxes are normalized to the 
maximum value reached during that particular time-interval. 
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Fig. 7. Left panels: Spectral patterns in the a-flux plane showing clockwise and counterclockwise trends. Right panels: Swift- 
XRT light curves at three different bands: soft (0.2 - 3.0 keV, red), medium (3.0-5.0keV, green), and hard (5.0 - 10.0 keV, red) . 
Capital letters are referenced in the discussion in the text. In both left and right panels, the fluxes are normalized to the maximum 
value reached during that particular time interval. 



Flare 4) This flare has three patterns: CCW,CW, and CCW. It 
starts with a CCW loop. Initially there is a spectral change 
with an almost steady flux (A), followed by a weak spectral 
hardening with a flux increase (B). The CW loop is domi- 
nated by the cooling time. The final CCW pattern is char- 
acterized by a flux increase and a very rapid spectral hard- 
ening. Also in this case the new flaring component seems 
to start from the hard band. The flux decrease (C) is almost 
achromatic in this case, and the escape time may be domi- 
nant over the cooling one ( T esc <K r coo ;). 

Flare 5) This flare has three patterns: CW,CCW, and CW. It 
does not have features different form the previous ones. 

CW loops were ob served for the same sourc e in the past by 



CW loops were ob served tor the same source in the past by 
Takahashi et~ai1dl996l) with ASCA observations. lZhangl (l2002) 



observed CCW in BeppoSAX observations on 21 April 1998. 
Simi lar kind of loops w ere also observed in PK S2 155-304 
bv iKataoka et al.1 d2000h and 1H1426+428 by iFalcone et af 
d2004 . 



6. Spectral parameter trends. 

We compare the trends among the spectral parameters with 



those re sulting from the statistical analysis in ITramacere et al 



d2007bl) . We extend the data set presented in that work with the 
results from the observations analysed in this paper. It is worth 
noting that the whole data set spans about ten years, sampling 
the source in very different brightness states, making the final 
result more statistically significant. 



6.1. E p -b trend 

The evolution of the curvature parameters as a function of the 
peak energy points out relevant features of the acceleration 
process. This analysis indicates that as the peak energy of 
the emission increases, the cooling time scale shortens and 



can compete with the time scales for acceleration. In this 
case, it is possible to observe a bias in the E p - b relation, 
due to the cooling time scale dominating over time scales for 
acceleration. 

We analysed the data in Tab. [3] identifying cooling-dominated 
observations (lines labeled with (c)), characterized by a strong 
flux decrease and strong spectral softening (a =* 2). In Fig. [8] 
we plot with green empty circles the whole XRT data set from 
Tab. [3] and with black squares the points without the strong 
cooling contamination discussed above. 



In the same figure we report also data from ITramacere et al 
(2007b), showing that the trend in our sample is consistent 
with that from the historical data. This agreement confirms that 
the spectral curvature decreases as E p moves toward higher 
energies. 

This phenomenon can be interpreted according to two dif- 
ferent scenarios. A first scenario is that of an energy-dependent 
acceleration p robability process (EDAP). Within this context, 
Massa ro et al. d2004 showed that for acceleration efficiency 
inversely proportional to the energy itself, the energy distri- 
bution approaches a log-parabolic shape. According to this 
model, the curvature (r) is related to the fractional acceleration 
gain (e) by r oc — !— . A possible example is given by particles 

lOg £ 

confined by a magnetic field, whose confinement efficiency 
(Pace) decreases as the gyration radius (R L ) increases. From 
E p oc s and r oc — the negative trend between E p and b 
follows. This trend is in agreement with the observed trend. 



An alternative explanation is provided by the stochastic ac- 
celeration framework (SA), with the presence of a momentum- 
diffusion term. In this scenario, the diffusion term acts on the 
electro n spectral shape b roadening the distribution. In par- 
ticular, iKardashevI (I1962I) showed that a log-parabolic spec- 
trum results from a Fokker-Planck equation with a momentum- 
diffusion term and a mono-energetic or quasi -mono energetic 
particle injection. The results presented in ITramacere et al 
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O Swift Xrt Apr- Jul 2006 

■ Swift Xrt Apr-Jill 2006. nnn-cool iny dominated 
• Tramacere et al. 2007b (data from 1997 to 2005) 




Fig. 8. Scatter pl ot of the curvature (b) vs . E p . Red circles rep- 
resent data from lTramacere et al.l (I2007bl) spanning from 1997 
to 2005, instruments: ASCA, BeppoSAX, XMM-Newton. Black 
boxes represent Swift data from the present analysis without 
the cooling-dominated events. Empty green circles represent 
the whole XRT data set presented in this paper. 



Swift Xrt Apr-Jul 2006 

Tramacere et al. 2007b (data from 1997 to 20051 



1 E n (keV) 



Fig. 9. Scatter plot of L„ vs E p . Red circles represent data from 
Tramacere et al. (2007b) spanning from 1997 to 2005, instru- 
ments: ASCA, BeppoSAX, XMM-Newton. Blue boxes represent 
Swift data from the present analysis. The dashed dotted line 
represents the best fit using a power-law and the dashed line 
represents the best fit using a broken power law. 



2007b) rely on the theoretical prediction from the iKardashevI 
19621) model, that the curvature term (r) is inversely propor- 



tional to the diffusion term (D): 



1 



Dt 



(4) 



This relation leads to the following trend among the peak en- 
ergy of the electron distribution (y p ), the synchrotron curvature 
(b = r/5), r, and E„: 



\n(E p ) = 2ln( 7p ) + 3/(5b). 



(5) 



In paper II ((Tramacere 2009) we will discuss in detail how 
the stochastic acceleration can be used to reproduce this trend 
drawing a physical scenario that fits this phenomenological pic- 
ture. Here, we just remark that both the momentum-diffusion 
term D and, the fractional acceleration gain e can explain the 
anticorrelation observed between E p and b. 

6.2. Sp(L p )-E p trend 

The trend between S p and E p provides interesting indications 
concerning the driver of the spectral changes in the X-rays, in 
terms of the synchrotron emission. This trend can be used to 
understand how the luminosity of the jet evolves as the particle 
energy increases. In the fra mework of the synchrotron theory 
(Rvb icki & Lightmanll 19791) . the dependence of S p on E p can 
be expressed in the form of a power-law: 



S p *E a p . 



(6) 



In fact, starting from the following functional relation the 
SED peak height reads: 



S P <* n(y 3p ) y 3 



B 2 6 4 



(7) 



and the peak energy is given by: 



e p 00 Ji P B 6 



(8) 



where yi p is the peak of n(y) y 3 , B the magnetic field, and 
6 is the beaming factor. 

If we take into account a log-parabolic distribution for the elec- 
trons emitting the X-ray photons, then it is easy to show that 
the total emitters number N - J n(y)dy n{y p ) y p , with y p 
the peak of n(y) y. This in terms of S p implies 



ocNylB 1 5 A . 



(9) 



Thus a = 1 .0 applies when the spectral changes are dom- 
inated only by variations of the electron average energy and 
N is constant, a = 1.5 applies when the spectral changes are 
dominated by variations of the average electron energy but N 
is not constant; a — 2 as for changes of the magnetic field; 
or = 4 if changes in the beaming factor dominate; formally, 
a = oo applies for changes only in the number of emitting par- 
ticles, which implies either a change in the electron density or 
a change in the source size. The case a = 1.0 holds only if 
the electron distribution is a log-parabola and if the curvature 
is constant. 

In order to have a deeper understanding of the relation to the 
jet energetics, we plot on the y axis of Fig. [9] L p — S p An D 2 L , 
where Di 134.1 Mpc is the luminosity distanc^H We report 
the L p values both in Tab. [2] and [3] The resulting scatter-plot 
represents the L p - E „ trend obtained from m erging our data 
sample with that from lTramacere et al.l d2007bh . 
We fitted the data by means of a simple power law (PL) L p oc 
E^, and a broken power-law (BPL) 



1 We used a flat cosmology model with Ho = 0.71 km/(s Mpc) 
n M = 0.27 and Q A = 0.73. 
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L p oc E" 1 , E p < Eb 



(10) 



E„ > Eb. 



The simple power-law fit gives a value of a — 0.42 + 0.06, 
and the broken power law gives a break energy Eb = 1 + 1 keV, 
with the two slopes a\ = 1.1 ± 0.2 and a 2 = 0.27 + 0.07. 
We focus on the results from the BPL fit. Although the break 
energy is not well constrained, it is worth to note that spec- 
tral slopes are quite different with a high statistical signifi- 
cance. This break in the trend implies that for E p < IkeV and 
L p < 10 45 erg/s the driver follows the relation with a 1.0 (we 
define this state the quiescent sate), whilst for E p > IkeV and 
L p > 10 45 erg/s, the driver is ruled by a 0.2 (we define this 
as the high state). 

In Paper II dTramacerd l2009) we will discuss in detail how the 
stochastic acceleration or the energetics of the jet can be used 
to reproduce this trend. 
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Fig. 10. Three different spectral shapes from the data presented 
in this paper. Red boxes represent a power-law spectrum ob- 
served during Orbit 3 (ObsID 00030352010) on 06/16/2006. 
Blue diamonds represent a spectrum that is a log-parabola 
with a low energy power-law tail, from Orbit 3 (ObsID 
00030352005) on 04/25/2006. Green circles represent a log- 
parabolic spectrum from Orbit 4 (ObsID 00030352008) on 
06/14/2006. 



7. The Swift-UVOT Swift-XRT connection: the 
low-energy power-law tail in the electron 
distribution. 

Both E p - b and S p - E p trends allowed us to understand the 
shape of the electron distribution for particles emitting at ener- 
gies close to E p . Anyway, the particle energy distribution can 
develop a low energy shape quite different from the extrapola- 
tion of the high energy branch. This difference is relevant for 
discriminating among different acceleration processes. In this 
perspective, the connection between the UV and X-ray spec- 
tra can give useful information about the low energy tail of the 



electrons emitting in the X-ray. We analysed carefully all of the 
spectra with simultaneous X-ray and UV observations. As re- 
sult we found that joint UVOT-XRT SEDs can be classified in 
three categories: 

a) described by a log-parabola (LP) 

b) described by a power law (PL) 

c) described by a spectral law that is power law in the low 
energy tail, t urning into a log-par abola in the high energy 
one (LPPL) ( Massaro et al. 20061) . whose functional form 
can be expressed by 



vF(v) = N(v/v c T a \ v<v c 
vF(v) = N(v/v c y (a ^ b l0S{vM \v>v c 



(11) 



where a v is the spectral index of the SED (v F (v)), and v c 
is the frequency where happens the turn-over in the SED ( 
Fig. El. 



From the analysis of this spectral behaviour it is possible to 
constrain the minimum energy of the radiating electrons. In fact 
electrons radiating mainly in the UV band have a Lorentz factor 
juv satisfying the following condition dRvbicki & Lightman 
19791) : 



10 15 ffz « 3.7 xlO 6 B6y uv 2 /(l+z) 
1.6 x 10 4 . 



(12) 



yuv 



1+z 



B6 



If the spectral shape is consistent with the same log-parabola 
extending from the X-ray band down to the UV band (case a), 
then it means that electrons radiating at UV frequencies belong 
to the same electron population and according to Eq. [12] we 
have y min < y uv . 

The condition y m - m > juv may occur when we observe a 
PL (case b) or a LPPL (case c). In fact, if y min > y uv then 
the spectra in the X-ray-to-UV band will be described by the 
asymptotic low-energy approximation of the single particle 
synchrotron e mission that is a power law w ith slope a v -4/3 
(SED oc v 4/3 ) dRvbicki & Lightmanlll979l) . 
For both case b) and c) we note that our data give a v ^ 
[0.25 - 0.4], a value very different from the asymptotic syn- 
chrotron kernel expectation. This hints that both cases the UV 
photons are likely to originate from an electron distribution 
that has a power law tail in the energetic range radiating in 
the UV-to-soft-X-ray band. A phenomenological option to ex- 
plain the case c) is an electron distribution that is a power law at 
low energies with a log-parabolic high-energy branch (LPPL) 
(Ma ssaro etal.ll2006l) : 



«(r) = K (y/y c ) s , y <y c 

n(y) = K{yly c )- {s+rLog{yiy ' ) \y>y c 



(13) 



where y c is the turn-over energy. 

In the case b), the electron distribution is assumed to be a pure 
power-law. 

Interestingly for case b) and case c) we can also constrain 
the typical slope of the power-law branch of the electron dis- 
tribution, using the well known relation between the spec- 
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tral index in the particle d istribution s and that in the SED 
dRvbicki & Lightmanll979T) : 



SED oc v 



„-(»-3)/2 



(14) 



For the typical values of a v observed in our data set, the 
resulting value of s is in the range s =s [2.2 - 2.5]. 

The presence of a power-law feature and the range of ob- 
served spectral indices are relevant both in the context of Fermi 
first-order acceleration models and from an observational point 
of view. 



for this spectral shape is a power-law at low energy with a log- 
parabolic high-energy branch (Eq. [Oi l. On the contrary, the 
03/31/2005 SED can be modeled using a log -parabolic elec- 
tron distribution that we express in terms of the peak energy as 
(LPEP): 



n{y) = K 10-'' (logWy " ))2 . 



(15) 



From the obse rvational side , it is worth to remark that Waxman 
( 1997 ) and Meszaros| d2002l) . studying the the afterglow X-ray 
emission of y-ray bursts (GRB), inferred an electron distribu- 
tion index s 2.3 +0.1. This is very close to those found in 
our data, but coming from a quite different class of sources. 

From a theoretical point of view, there are several 
works concerning relativistic shock acceleration models 
that start from different analytical o r numerical approaches 
and find values of s * \2.2 - 2.4] dAchterberg et alJbOOlt 



Since we do not know the actual shape and flux at TeV 
energies (we have only an estimate of the flux at ^ 1 TeV) we 
cannot constrain in d etail the SSC model u sing the canonical 
B - 5 plane analysis dTavecchio et al.ll 1998D . Yet, we can still 
use the results reported in Sect. 3 and 7 to constrain some of 
the SSC parameters. 

The first constraint comes from the source variability (see Sect. 
3): 



R/6< 1 x 10 14 cm. 



(16) 



Gallant et al.l 119991 ILemoine & Pelletierl 120031 

2005j[ lEllison & Double! |2004 . These values are consistent y„„„ < 1.6 x 10 4 
with those from our data set. 



A second constraint comes from the analysis of the UVOT con- 
nection with the XRT data that gives an upper limit to y m ,„ 



Blasi & Vietri 



l+z 
B5 ' 



(17) 



The power-law feature is also consistent with a purely stochas- 
tic scenario. Anyway, the usual limitation of the stochastic 
model to explain a universal index relies on the fine tuning 
required on the ratio of the acceleration time scale to the loss 
time (s 1 + tacc/tesc), in order to match the observed values. 
Moreover, we stress that a power-law electron distribution 
n(y) oc y~ 2 - 3 is not compatible with a Maxwellian-like dis- 
tribution (n(y) oc y 1 ) resulting from the equilibrium of SA 
processes without relevant particle escape. 

In conclusion both case c) and b) are better explained 
better by a first order process. We will discuss this topic further 
in Sec. 9. 



8. SED modeling and GeV/MeV predictions. 

We model the SEDs of the three observations that have simul- 
taneous XRT, BAT and UVOT data, using a standard one-zone 
SSC scenario. The only useful TeV data found in the litera- 
ture are from Wh ipple observations on June 18,19, and 21 
(Lich ti et al . 2008). These data, almost simultaneous only with 
the 06/23/2006 Swift pointing, provide only the TeV flux, with- 
out giving a description of the spectrum. For this reason they 
are used only to estimate the TeV flux level during that point- 
ing. 

To have a feeling of the spectral and flux range of 
variability, we also plo t Swift-XRT data from 03/31/2005 



To constrain the value of the maximum electron energy we can 
use the maximum energy of the synchrotron emission, taking as 
the corresponding energy the most energetic bin of the BAT de- 
tector with a significant signal (^ 50 keV) 



y max > 1.8x10° 



l+z 
B6 



(18) 



A further constraint on y max can be set by estimating the 
typical electron energy req uired to produce TeV photons 
(Bed narek & Protherodl 19991) . Based on the historical data for 
Mrk 421, the spectrum typically reaches up to 10 TeV, and this 
implies 



Jmax > 2 X 10 



7 1 +Z 



(19) 



And following iBednarek & Protheroe (1999), combining Eq. 
[l9]and[T8]we obtain an upper limit on the magnetic field: 



B/6 < 0.008 G. 



(20) 



(Tra macere et al.l l2007al) . and some TeV SED s repre s enting 



the source in d ifferent flaring states dAlbert et al. 2007 
Yadav et al.ll2007h (see left panel of Fig. E) 



The 2006 SEDs that we want to model have a power-law 
spectral dependence between the UVOT and XRT bands. As 
described in Sect. 7, the most generic distribution accounting 



The best fit model was obtai ned by combining our numer- 
ical SSC code dTramacerell2007 ) with a numerical minimizer. 
According to the observationally derived constraints, we fixed 
the value of the beaming factor at 6 = 25, the magnetic field 
B - 0.1 G, the source size R - 2.1 x 10 15 cm, and tuned 
Jmin - 1100 to get the right Compton dominance, leaving 
r, s, y c , Jmax and N as free parameters. The resulting best fit 
parameters are reported in Tab. [4] and [5] 

In the left panels of Fig. QT| we show the best fit results 
for the SSC model for the 2006 data and for the 31/03/2005 
pointing. In the right panels we show the corresponding elec- 
tron distributions. The values of the electron curvature are 
consistent with those observed in the X-ray emission accord- 
ing to the relation b r/5. The typical value of the ratio 




Fig. 11. SSC fits of three different observations with simultaneous UVOT XRT and BAT data. Left panels show the SSC model, 
from top to bottom: solid circles represent data from 04/22/2006, 06/23 /2006, 07/15/2006 Swift observations. Green triangles 
show Swift-XRT data on 03/31/2005, data from lTramacere et al.l d2007bl) . Solid grey polygons represent no n-simultaneous EBL 
corrected TeV data. Solid gray squares represent the high state on 2001 observed by Whi pple, data from | Albert et alj (I2007I) . 
Solid gray diamonds represent the average 2004-2005 TeV spectrum as observed by MAGIC dAlbert et al.l200 7 |). Solid gray righ t 
triangles represent the average spectrum from December 2005 to February 20 06 as observ e d by TACTIC (lYadav et al.ll2007l) . 
The solid red triangle represents a Whipple observation on June 18,19,21 from lLichti et al.l (120081) . that is very close in time to 
our 06/23/2006 data set. The solid lines represent the best fit by a SSC model to our simultaneous Swift observation, and the 
dashed line is the best SSC fit to the Swift data on 03/31/2005. The dotted lines represent the modeling of the galaxy contribution 
by means of black body spectral shape. Right panels show the electron distributions for the SSC models in the left panels. The 
solid lines in the right panels represent the electron distributions for the best fit models of the three 2006 Swift observations, the 
dotted lines represent the extrapolation of the LP branch of the LPPL distribution, and the dashed lines represent the electron 
distribution for the 03/31/2005 data. 



Ue/uB — [120 - 180] results in a particle dominated jet in 
agreement with the typical values for TeV HBLs peaking in 
the hard X-ray dKino et al.ll2002t ISato et alj|2008h . The low- 
energy power law branch of the electron energy distribution 
follows n{y) oc y -23 , in agreement with the analysis presented 
in Sect. 7. 

We remark that UV data are weakly contaminated by the 
host galaxy emission. 

Regarding this, we plot in the left panel of Fig. [TT] the 
contribution of the galaxy emission modeled as a black body 
spectral shape (dotted-line). This plot clearly shows how the 



UVOT data taken into account lie beyond the galaxy emission 
spectral cut-off. On the contrary, during the low states, the pos- 
sible power-law tail would lay typically at optical and lower 
frequencies, where the emission is typically dominated by the 
galaxy. This implies that we actually don't know whether the 
electron distribution has such a power-law tail also during the 
lower state. The SED model of data on 31/03/2005 shows in 
fact that the UVOT points are compatible with the synchrotron 
emission from a log-parabolic electron population. 
Moreover, a power-law low-energy branch if present in the 
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2005 SED would yield an index harder than those from 2006, 
but with much lower flux. 

Another interesting analysis can be performed by looking at 
the acceleration time scales determined by the rate of the vari- 
ation of Ep. In particular, we can use the e-folding time of E p 
to infer the electron acceleration time scale r acc . In detail, from 
E oc y 2 it follows that the rate of change of E p can be linked to 
the rate of change of y for electrons radiating in the hard X-ray 
band. In Fig [12] we show a flare having well constrained val- 
ues of E p for the rising side. In this figure we plot the values 
of E p as a function of time. The e-folding time for E p results 
in a y e-folding time of about 0.15 days, which is the value of 
the typical acceleration time scale. This is the time scale that 
will compete with the cooling one to generate the equilibrium 
in the n(y) distribution. In particular taking into account the 
synchrotron cooling: 



t coo1 - 107(5%) s 



(21) 



the equilibrium of the distribution for a log-parabolic or 
Maxwellian-like distribution is the peak of n{y). In the case 
of LPPL distribution it will be constrained between y c and y p p . 
Imposing t coo1 = T acc we find: 



Jeq = 



10 8 



B 2 5 



(22) 



For T acc — 0.15 days, 6 — 25 and using y eq of the order of 10 4 
the resulting value of the magnetic field is 



B 



10 8 



Tacc&Jp 



0.16 G 



(23) 



which is consistent with the upper limit of Eq. [20]and with our 
best fit field strength. 
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Fig. 12. Lower panel: the increase of E p during a flare starting 
on 06/14/2006, fitted by an exponential function. Upper panel: 
the corresponding HR. 



8.1. The connection between the UV slope and that at 
Me V/Ge V energies. 

A further and very relevant consequence of the spectral shape 
in the UV band is its connection with the shape expected at 
MeV/GeV energies, and in particular in the Fermi Gamma-ray 
Space Telescope-LAr band. 

In order to explore this we simulated SSC emission accounting 
for a typical HBL SED. We used as the electron distribution 
a power-law at low energy with a log-parabolic high energy 
branch (LPPL) (Eq. [13) . The SSC parameter are reported in 
the caption of Fig. [13] In the left panel of this figure we show 
different SEDs obtained by varying the s parameter in the range 
[2.1:2.5], and in the right panel we show how the spectral shape 
in the UVOT band is tightly correlated with that in the Fermi- 
LAT energy range. This means that a comparison between the 
two spectral slopes performed with simultaneous observation 
can hint whether actually the UV photons from the same com- 
ponent emitting in the X-ray are responsible for the SSC emis- 
sion. 



9. Discussion. 

In this paper we presented analyses of the spectral and 
flux evolution of Mrk 421 during the spring/summer 2006 
Swift observations, with the June pointings monitoring the 
source almost continuously for 12 days. 

During this period the source exhibited both flux levels and 
SED peak energies at their historic maximum until 2006. 
This very intense flaring state of Mrk 421 represented a 
unique opportunity to test the correlations amo ng the spectral 



Pi 
2( 



amme ters presented in previous works dTramacere et al 
2007bl) . expanding the temporal spanning to about 10 years 
and enlarging the volume of the parameter space. 



The spectral evolution of the source and the patterns 
in the a-flux plane (Sect. 4) suggest that each flare has a 
different characterization in terms of competition among the 
relevant time scales. Some flares showed a rise in the hard and 
medium X-ray band much faster than that observed in the soft 
X-rays. This behavior is likely to be explained with the flaring 
component starting in the hard X-ray band and may suggest 
that the driver of those flares is the rapid injection of very 
energetic particles rather than a gradual acceleration. 



The S 



p 



Tramacere et al 



and E 



(2007 



- b trends follow those presented in 
with E p and S p correlated and E p ,b 
anti-correlated. These trends are relevant in order to understand 
the physical mechanisms driving the evolution of the electron 
distribution under the effect of acceleration proc esses and may 



represent a common scenario for HBLs. In fact. iMassaro et al 



(2008) confirmed that the E p - b relation holds for five of 



the TeV HBLs included in their analysis (PKS 0548-322, 
1H 1426+418, Mrk 501, 1ES 1959+650, PKS 2155-304)), and 
that as for the S p — E p the only exception to the previous list is 
PKS 2155-304 

The E p - b relation shows a possible signature of acceler- 
ation processes leading to curved electron distributions, with 



A. Tramacere,ef al: Swift observations of the very intense flaring activity of Mrk 421 during 2006. 



15 















































\ 


~ UVOTband 


l\\ w 


lOGeV band||j 



- 

8,0.35 



Log(v) (Hz) 



0.25 - 
0.25 



0.35 0> 
UVOT spectral slope 



Fig. 13. The strong correlation between the UVOT spectral slope and that in the LAT band for a SSC scenario. In the left panel 
we show the SSC SEDs and in the right panel we show the correlation between the spectral slope in the UVOT band and 
that in the Fermi-LAr band. Basic parameters for the SSC model are: particle number density N — 10 cmT^ \ magnet ic field 
intensity B = 0.1 G, minimum electron energy y m - m = 10 3 , maximum electron energy y„ mx — 5 x 10 6 , emitting region radius 
R = 3.0 x 15 cm, beaming factor S = 15, electron curvature r — 1.0, power law of index of the low-energy tail of the electron 
distribution s = [2.2 : 2.5], and turn-over energy in the electron distribution y c — 5.0 x 10 4 

Table 4. SSC best fit results for the 2006 Swift observations and using as electron distribution a log-parabola with a power-law 
low-energy branch (Eq.fl3l, 



Date 


B 


R 






6 


N(*) 


r 


s 


y c (**) 




Yfi* 


**) 


Yma 


X 




ymin 




Ue/Ub 




G 


cm 








cm~^ 




























22-04-2006 


0.1 


2.1 


X 


10 15 


25 


13.5 


0.75 


2.30 


1.75 x 


10 s 


2.0 x 


10 3 


2.5 


X 


10 6 


l.l 


X 


10 J 


119 


23-06-2006 


0.1 


2.1 


x 


10 15 


25 


15.0 


0.65 


2.30 


2.85 x 


10 s 


4.8 x 


10 s 


4.0 


X 


10 6 


l.l 


X 


10' 


135 


15-07-2006 


0.1 


2.1 


X 


10 15 


25 


21.0 


0.85 


2.32 


2.50 x 


10 5 


1.0 x 


10 4 


3.0 


X 


10 6 


l.l 


X 


10 3 


185 



Table 5. SSC best fit results for 31/03/2005 Swift observation and using a log-parabolic electron distribution as defined in Eq. 

El 



Date 


B 


R 


6 


JVC*) 


r y P (**) 


/max 


Jmin 


Ue/Ub 




G 


cm 




cm" 3 










31-03-2005 


0.075 


1.5 x 10 15 


25 


4 


1.3 4.0 xlO 3 


2.5 x 10 6 


1.1 x 10 3 


220 



(*) N = J n{y)dy. 

(**) Do not confuse y c (the turn-over energy in Eq.|13t with y p (the peak energy in Eq.|15l). 

(***) Tp represents y p obtained from extrapolating the log-parabolic branch of the distribution described by Eq. [T3] These values can be 
compared with y p from 31/03/2005. 



the curvature decreasing as the acceleration gets more efficient. 
A first scenario supporting this picture is that of an energy 
depe ndent acceleration probability process ( iMassaro et al.l 
120041) . An alternative explanation is provided by the stochastic 
acceleration fra mework, with the p resence of a momentum - 
diffusion term (Kardashev 1962 : Tramacere et al. 2007b). 



Both scenarios predict a negative correlation between E p and 
b and are consistent with the the data plotted in Fig. [8] It is 
worth noting that the data presented here substantially stretch 
the parameter space, confirming that the relation between the 
curvature and the peak energy follows the same trend from 
the faintest to the strongest flaring activity, hinting that the 



acceleration process is actually be the same. 
A future investigation may be focused on understanding which 
are the physical parameters that can tune the acceleration 
process. In this regard, an interesting investigation may be the 
connection of the stochastic scenario with the turbulence spec- 
trum of the MHD waves and the resulting d i ffusion coefficient 
(|Park & Petrosianll 19951: iBecker et alj|2006l : iKatarzvnski et al 



2006|; iStawarz & Petrosianl 120081) . Understanding the role 



of the turbulence spectrum, and in general, the role of the 
fluctuation on the particle energy gain, is a complex task and 
requires a more detailed analysis that is beyond the purpose of 



16 



A. Tramacere,ef al: Swift observations of the very intense flaring activity of Mrk 421 during 2006. 



this paper and will be presented in Paper II ( Tramacerell2009 ). 



The S p - E p trend demonstrates the connection between 
the average energy of the particle distribution and the power 
output of the source. The expected power-law dependence 
S p oc E a p compared to the observed values of a, rules out 
B, 6 and N, indicating y p as the main driver of the S „ — E p 
trend, confirming the result from iTramacere et al.1 d2007bl) . A 
more detailed analysis of the scatter plot reported in Fig. [9] 
revealed that the trend has a break at about 1 keV, where the 
typical source luminosity is about L p 10 45 erg/s. This break 
may be interpreted as the marker of the competition between 
systematic and momentum-diffusion acceleration or in terms 
of energetic content of the jet. We will study such a scenario in 
Paper II, where we will take into account also the role E p - b 
trend on S p - E p . 

Another interesting result from our analyses is the presence 
a power-law tail connecting the UV data to the soft-X-ray. As 
explained in Sect. 7, we cannot determine whether this tail is 
a characteristic only of highly energetic flares, or whether it is 
present during low states too. The puzzling aspect in the latter 
case would be to reconcile a harder index with a lower flux. 
Anyway, if this feature is present only during strong flares 
then it may be argued that a relevant change has taken place 
in the acceleration environment. In fact, both SA and EDAP 
scenarios require the presence of a significant escape term to 
develop a power-law tail. 

The value of the power-law spectral index in the electron 
distribution s 2.3 is very close to the pre diction from rela- 



tivistic Fermi first ord er accele ra tion mod els (Achterberg et al 



2001; Galla nt et al l 1 19991 iLemoine & Pelletierl 120031: 



Blasi & Vietril 120051: lEllison & Double! l2004 and it is 
similar to that found in the X-ray observations of GRB 



afterglows ( Waxmanlll997 ). 

This observational feature, which is in nice agreement 
with first order acceleration models, may suggest a kind of 
contradiction with the results from the E p - b trend, supporting 
a stochastic scenario. A more careful analysis shows that 
such a cont radiction is only apparent. In fact, as pointed out 
recently by ISpitkovskvl (12008b . diffusive shock acceleration 
models rely on the presence of a magnetic turbulence near the 
shock responsible for the scattering process, but these models 
do not clearly determine the possible role of the turbulence 
in the acceleration. The indication from our observational 
analysis may be simultaneous roles for the first and second 
order processes, with the stochastic acceleration arising from 
the magnetic turbulence and signed by the curvature and more 
generally by the E p - b relation, and the first order signed by 
the slope of the electron power-law tail. 

In such a circumstance the study of the curvature may offer an 
observational constraint for the turbulence properties, and may 
be recursively used to constraint the scattering process near the 
shock in the first order models. 



our phenomenological picture. The result from this numerical 
study is a relativistic Maxwellian, and a high-energy tail with 
s = 2.4 + 0.1, plus an exponential cut-off moving to higher 
energies with time of the simulation. 

This phenomenological scenario underlines the relevance 
of the multi-wavelength observations, and warns about the 
pitfalls of extrapolating the observed spectral shape over too 
wide a range. In this regard we stress the unique capabilities of 
Swift to perform simultaneous UV-to-X-ray observations. 

A further point concerning the UV observations is given by 
the expected correlation of this spectral band with MeV/GeV 
band. In fact due to the effect of the Klein-Nishina suppression 
in the Inverse Compton process, we expect the UV photons 
to be the most efficiently up-scattered at GeV energies by 
the electrons radiating in the X-ray. This makes a strong 
correlation between the spectral slope in the UV-to-soft-X-ray 
band with the slope at MeV/GeV energies and in particular in 
the Fermi-LAr band. A test of such a correlation would be 
useful in order to understand whether the photons up-scattered 
at y-ray energies are co-spatial with the electrons emitting in 
the X-ray. 



10. Conclusions. 

In conclusion the present work shows the complexity of the 
physical scenario at work in the jets of HBL objects. We have 
shown that the flaring activity of Mrk 421 not only causes 
huge flux variation but also results in complex spectral evo- 
lutions and drastic changes in the electron energy distribution. 
The evolution of the spectral parameters, in particular E p and 
b, agrees with an acceleration model where the curvature term 



is inversely proportional to E p , which is consistent with both 
SA and EDAP scenarios. This supports the hypothesis that the 
spectral curvature is related to the acceleration rather than to 
the cooling process. In fact we find less curved spectra as E p 
increases, leading to the conclusion that the cooling is more 
relevant to determine the reaching of the equilibrium energy 
rather than to determine the spectral bending. The presence of 
a power-law low-energy tail remains puzzling because we can- 
not determine whether it develops only during strong flares or 
whether it is a common feature. The value s 2.3 is very close 
to the universal index from relativistic shock models and could 
be used to understand the relative weight of the first order pro- 
cess during strong flares, relative to the diffusive process in- 
voked to explain the spectral curvature. We remand these open 
questions to Paper II, were we will focus mainly on the theoret- 
ical interpretation of the phenomenological picture presented 
here. 

As a final remark, we stress that as already pointed out in 
Mas saro et al. ( 2006h . the curvature observed in the X-ray spec- 



We n ote moreover that the model presented by ISpitkovsky 
(2008), in which the shock acceleration process is studied in a 
self-consistent fashion, gives results that are compatible with 



tra of HBLs is reflected in the Te V spectrum. The consequence 
is that the TeV cut-off observed in the nearby HBLs is al- 
most entirely due to the intrinsic electron curvature rather than 
to the interaction with EBL photons. This scenario is there- 
fore consistent with a low EBL density that makes the uni- 
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verse more transparent to high-energy radiation than previously 
assumed, in agreement with the discovery of HBLs objects 
at large re dshift like H 2356-30 9 (z=0.165), 1ES 1101-232 
(z=0.186) jAharonian et ai] |2006l) and 1ES 128-304 (z=0.182) 
dAlbert et al.ll2006l) . 
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Table 2. Swift-XRT Spectral analysis of Mrk 421 . 



Orbit a 


b K E p * E p 
keV keV 




flux 0.3-10 keV L p 

10 ^erg cm~~ s - ' 10 45 erg/s 


tf/dof 




Obsld 00206476000 


Date 4/22/2006 MJD 53847.252792 









01 


1.69(0.02) 


0.09(0.05) 


0.347(0.005) 


> 10 


- -(- -) 


- -(- ") 


2175.7 


- -(- -) 


0.861(137) 


02 


1.66(0.02) 


0.13(0.04) 


0.349(0.004) 


> 10 


- -(- -) 


- -(- -) 


2191.2 


- -(- -) 


1.181(197) 


03 


1.71(0.02) 


0.07(0.04) 


0.336(0.004) 


> 100 


- -(- -) 


- -(- -) 


2119.3 


- -(- -) 


1.211(196) 


04 


1.67(0.02) 


0.12(0.04) 


0.334(0.004) 


> 10 


- -(- -) 


- -(- -) 


2093.5 


- -(- -) 


1.224(170) 


05 


1.67(0.02) 


0.12(0.04) 


0.369(0.004) 


> 10 


- -(- -) 


- -(- -) 


2299.0 


- -(- -) 


0.902(207) 


06 


1.7(0.02) 


0.17(0.04) 


0.336(0.004) 


7(5) 




727(34) 


1980.1 


1.56(0.07) 


1.049(192) 


07 


1.69(0.03) 


0.21(0.06) 


0.306(0.005) 


5(4) 




635(32) 


1761.1 


1.37(0.07) 


1.016(94) 


08 


1.86(0.03) 


0.15(0.07) 


0.275(0.006) 


3(2) 


2 9 +3 

-0.8 


475(16) 


1452.2 


1.02(0.03) 


1.025(58) 


09 


1.84(0.03) 


0.14(0.06) 


0.25(0.004) 


4(3) 


3 6 +3 


444(14) 


1347.2 


0.96(0.03) 


1.116(96) 


10 


1.79(0.03) 


0.15(0.06) 


0.247(0.004) 


5(5) 


5 +6 

—2 


465(22) 


1367.5 


1(0.05) 


0.935(93) 


11 


1.82(0.03) 


0.08(0.07) 


0.241(0.005) 


> 10 


- -(- -) 


- -(- -) 


1374.0 


- -(- -) 


1.163(80) 


12 


1.8(0.03) 


0.08(0.06) 


0.253(0.004) 


> 10 


- -(- -) 


- -(- -) 


1447.7 


- -(- -) 


0.991(98) 


13 


1.78(0.02) 


0.09(0.05) 


0.269(0.004) 


> 10 


- -(- -) 


- -(- -) 


1562.1 


- -(- -) 


1.020(105) 


14 


1.79(0.02) 


0.08(0.05) 


0.267(0.004) 


> 10 


- -(- -) 


- -(- -) 


1558.9 


- -(- -) 


0.989(127) 


15 


1.83(0.02) 


0.06(0.05) 


0.253(0.004) 


> 10 


- -(- -) 


- -(- -) 


1445.0 


- -(- -) 


1.166(119) 


16 


1.81(0.02) 


0.07(0.05) 


0.269(0.004) 


> 10 


- -(- -) 


- -(- -) 


1557.7 


- -(- -) 


1.103(122) 


17 


1.67(0.03) 


0.2(0.06) 


0.285(0.005) 


7(7) 


7+8 

'-2 


630(47) 


1691.3 


1.4(0.1) 


0.892(90) 


18 


1.85(0.03) 


0.12(0.06) 


0.25(0.005) 


5(6) 


- -(- -) 


- -(- -) 


1365.0 


- -(- -) 


0.863(89) 


19 


1.81(0.03) 


0.14(0.06) 


0.246(0.005) 


5(6) 


- -(- -) 


- -(- -) 


1362.0 


- -(- -) 


1.218(85) 


20 


1.83(0.02) 


0.12(0.05) 


0.251(0.004) 


5(6) 


- -(- -) 


- -(- -) 


1382.7 


- -(- -) 


0.945(115) 


21 


1.81(0.02) 


0.17(0.05) 


0.276(0.004) 


4(2) 


3.8!?, 


502(14) 


1500.4 


1.08(0.03) 


1.083(117) 


22 


1.83(0.02) 


0.09(0.05) 


0.294(0.004) 


9(15) 


- -(- -) 


- -(- -) 


1649.9 


- -(- -) 


1.127(130) 


23 


1.77(0.02) 


0.12(0.05) 


0.275(0.004) 


9(13) 


- -(- -) 


- -(- -) 


1586.6 


- -(- -) 


1.044(122) 


24 


1.82(0.02) 


0.09(0.05) 


0.267(0.004) 


9(15) 


- -(- -) 


- -(- -) 


1505.8 


- -(- -) 


1.066(113) 


25 


1.8(0.02) 


0.12(0.05) 


0.283(0.004) 


7(8) 


7+12 

1 -2 


548(34) 


1589.9 


1.18(0.07) 


0.901(132) 


26 


1.8(0.02) 


0.12(0.05) 


0.282(0.004) 


6(7) 


- -(- -) 


- -(- -) 


1579.1 


- -(- -) 


0.994(125) 


27 


1.74(0.03) 


0.14(0.07) 


0.283(0.006) 


8(14) 


- -(- -) 


- -(- -) 


1636.3 


- -(- -) 


1.120(64) 



Obsld 00030352005 Date 4/25/2006 MJD 53850.267940 



01 1.84(0.01) 0.15(0.02) 0.273(0.002) 3(1) 3.3^ 481(5) 1468.4 1.03(0.01) 1.360(288) 

02 1.86(0.01) 0.16(0.02) 0.264(0.002) 2.7(0.7) 2.7+°* 454(4) 1391.7 0.977(0.009) 1.248(285) 

03 1.86(0.01) 0.15(0.02) 0.269(0.002) 2.9(0.8) 2.9+g;| 464(4) 1427.2 0.998(0.009) 1.090(287) 

Obsld 00030352006 Date 4/26/2006 MJD 53851.146539 

01(*) 1.86(0.02) 0.17(0.03) 0.42(0.004) 2.5(0.7) 2.5^ 715(9) 2192.5 1.54(0.02) 0.989(208) 

02 1.83(0.01) 0.14(0.03) 0.285(0.002) 4(2) 4.2+ 2 18 515(9) 1556.8 1.11(0.02) 1.257(253) 

03 1.82(0.01) 0.17(0.03) 0.294(0.002) 3(1) 3.4+°] 525(7) 1584.4 1.13(0.01) 1.380(254) 
Obsld 00030352007 Date 4/26/2006 MJD 53851.951963 



01 1.87(0.01) 0.16(0.03) 0.265(0.002) 2.6(0.6) 2.6+g;| 451(4) 1388.5 0.97(0.009) 1.140(258) 

Obsld 00030352008 Date 6/14/2006 MJD 53900.016100 



01 


2(0.03) 


0.12(0.07) 


0.25(0.005) 


1(0.2) 


i +02 
1 -0.4 


400(8) 


1227.0 


0.86(0.02) 


0.924(178) 


02 


2.07(0.03) 


0.14(0.06) 


0.236(0.004) 


0.6(0.2) 


6 +0 2 


384(10) 


1106.2 


0.83(0.02) 


1.036(173) 


03 


2.07(0.03) 


0.05(0.07) 


0.23(0.005) 


0.2(0.7) 


- -(- -) 


- -(- -) 


1130.2 


- -(- -) 


0.862(171) 


04 


2.01(0.03) 


0.2(0.07) 


0.236(0.005) 


0.9(0.1) 


9+0.1 

"• -0.2 


378(7) 


1111.7 


0.81(0.01) 


1.006(175) 


05 


2.1(0.03) 


0.2(0.06) 


0.213(0.004) 


0.6(0.2) 


6 +0 1 

" -0.2 


351(8) 


960.3 


0.76(0.02) 


1.089(175) 



(*) gti with biased exposure, the affected data are not taken into account in the analyses. 

The second third and forth columns report the best fit estimate for the model in Eq. [2] The fifth column reports the value of the SED peak 
analytically estimated from Eq. [2] according to the best fit results. The sixth and seventh columns report the E p and 5,, best fit estimates using 
as best fit model Eq.[3] In the eighth column we report the flux in the 0.3-10.0 keV band, evaluated by Xspec integrating the Eq.[2]model. The 
ninth colum reports the SED peak flux luminosity evaluated as L p - S p 4n D 2 L , where D L ^ 134.1 Mpc is the luminosity distance. In the last 
column we report the reduced^ 2 and the degrees of freedom concerning the Eq.[2]fit. 
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Orbit 


a 


b 


K 


E p * 
keV 


keV 


s P 


flux 0.3-10 keV 

10~ erg cm~ 2 


10 45 erg/s 


Xi/dof 


06 
07 


2.1(0.03) 
2.08(0.03) 


0.06(0.07) 
0.07(0.07) 


0.195(0.004) 
0.197(0.004) 


0.2(0.5) 
0.3(0.6) 


16 +0 1 
27 +0 2 

"••"-0.07 


343(18) 
332(14) 


939.1 
960.0 


0.74(0.04) 
0.71(0.03) 


1.099(147) 
0.873(143) 



08 2.05(0.04) -0.08(0.09) 0.192(0.005) -- --(--) --(--) 1028.5 --(--) 1.038(96) 

09 2.02(0.02) -0.01(0.04) 0.167(0.002) -- --(--) --(--) 872.3 --(--) 1.137(266) 



Obsld 00030352009 Date 6/15/2006 MJD 53901.490493 



01 


1.83(0.02) 


0.06(0.03) 


0.181(0.002) 


> 10 


- -(- -) 


- "(- -) 


1035.8 


- -(- -) 


0.979(270) 


02 


1.89(0.02) 


-0.02(0.04) 


0.162(0.002) 




- -(- -) 


- -(- -) 


941.5 


- -(- -) 


1.150(219) 


03 


1.87(0.03) 


0.01(0.06) 


0.14(0.002) 


> 100 


- -(- -) 


- "(- -) 


803.3 


- -(- -) 


0.842(155) 


04 


1.86(0.02) 


0.04(0.05) 


0.194(0.003) 


> 10 


- -(- -) 


- -(" -) 


1100.3 


- -(- -) 


0.936(200) 


05 


1.82(0.03) 


0.13(0.07) 


0.196(0.004) 


5(6) 


- -(- -) 


- -(" -) 


1078.0 


- -(- -) 


0.930(111) 


06 


1.88(0.03) 


0.11(0.07) 


0.191(0.004) 


4(5) 


- -(- -) 


" "(- -) 


1021.2 


- -(- -) 


1.313(97) 


07 


1.88(0.01) 


0.16(0.03) 


0.19(0.002) 


2.3(0.7) 


9 q+0.4 
il -'-0.3 


319(4) 


985.1 


0.686(0.009) 


1.017(273) 



Obsld 00030352010 Date 6/16/2006 MJD 53902.025296 



01 


1.92(0.01) 


0.08(0.03) 


0.201(0.002) 


3(2) 


3 2 +2 

D - -0.8 


337(5) 


1066.9 


0.73(0.01) 


1.287(222) 


02 


1.84(0.01) 


0.1(0.03) 


0.219(0.002) 


6(5) 


6!5 


404(13) 


1213.1 


0.87(0.03) 


1.228(237) 


03 


1.85(0.01) 


0.1(0.03) 


0.233(0.002) 


6(4) 




425(12) 


1283.5 


0.91(0.02) 


1.142(242) 


04 


1.83(0.01) 


0.07(0.03) 


0.254(0.003) 


> 


10 


- -(- -) 


- -(- -) 


1446.7 


- -(- -) 


1.234(226) 


05 


1.79(0.01) 


0.1(0.03) 


0.263(0.002) 


> 


10 


12 + ' 7 


546(32) 


1516.0 


1.17(0.07) 


1.146(263) 


06 


1.77(0.01) 


0.1(0.03) 


0.268(0.002) 


> 


10 


16+f 


590(43) 


1574.0 


1.27(0.09) 


1.227(263) 


07 


1.73(0.02) 


0.02(0.04) 


0.279(0.003) 


> 


100 




- -(- -) 


1773.0 




1.256(210) 


08 


1.76(0.02) 


-0.03(0.05) 


0.281(0.004) 








- -(- -) 


1808.0 




0.926(111) 


09(*) 


1.72(0.02) 


0.02(0.03) 


0.427(0.004) 


> 


100 




- -(- -) 


2741.8 




1.180(230) 


10 


1.69(0.02) 


0.04(0.03) 


0.261(0.003) 


> 


100 




- -(- -) 


1697.5 




1.110(220) 


11 


1.7(0.02) 


0.04(0.03) 


0.285(0.003) 


> 


100 




- -(- -) 


1833.8 




1.139(227) 


12 


1.73(0.01) 


0.03(0.03) 


0.302(0.003) 


> 


100 




- -(- -) 


1911.1 




1.258(240) 


13 


1.76(0.02) 


0(0.03) 


0.279(0.003) 








- -(- -) 


1764.0 




1.120(230) 


14 


1.74(0.01) 


0.04(0.02) 


0.29(0.002) 


> 


100 




- -(- -) 


1812.6 




1.034(288) 


15 


1.77(0.01) 


0.04(0.02) 


0.277(0.002) 


> 


100 




- -(- -) 


1677.1 




1.093(279) 


16 


1.81(0.01) 


0.05(0.03) 


0.246(0.002) 


> 


10 




- -(- -) 


1443.2 




1.128(242) 


17 


1.8(0.01) 


0.03(0.03) 


0.248(0.002) 


> 


100 




- -(- -) 


1485.6 




1.330(262) 


18 


1.84(0.01) 


-0.05(0.03) 


0.248(0.002) 








- -(- -) 


1512.6 




1.049(236) 


19 


1.84(0.01) 


-0.05(0.03) 


0.225(0.002) 








- -(- -) 


1383.1 




0.978(221) 


20 


1.81(0.01) 


-0.04(0.03) 


0.245(0.002) 








- -(- -) 


1517.3 




1.100(241) 


21 


1.7(0.02) 


0.03(0.03) 


0.287(0.003) 


> 


100 




- -(- -) 


1871.2 




1.038(240) 


22 


1.77(0.01) 


0(0.03) 


0.284(0.003) 








- -(- -) 


1769.6 




1.214(234) 


23 


1.75(0.02) 


0.03(0.03) 


0.293(0.003) 


> 


100 




- -(- -) 


1833.7 




1.116(226) 


24 


1.77(0.02) 


0.05(0.03) 


0.32(0.003) 


> 


100 




- -(- -) 


1930.8 




1.301(215) 



Obsld 00030352011 Date 6/18/2006 MJD 53904.038766 



01 


1.8(0.01) 


0.11(0.03) 


0.318(0.003) 


9(8) 


0+8 

-3 


633(27) 


1808.8 


1.36(0.06) 


1.036(286) 


02 


1.82(0.01) 


0.12(0.03) 


0.319(0.003) 


6(4) 


6!? 


600(15) 


1777.2 


1.29(0.03) 


1.227(286) 


03 


1.82(0.01) 


0.18(0.02) 


0.34(0.002) 


3.1(0.7) 


3-1^ 


601(6) 


1817.5 


1.29(0.01) 


1.343(317) 


04 


1.93(0.01) 


0.13(0.02) 


0.294(0.002) 


1.9(0.4) 


I n+0.2 


482(4) 


1511.3 


1.037(0.009) 


1.211(297) 


05 


1.94(0.01) 


0.11(0.03) 


0.292(0.002) 


1.9(0.5) 


1 9+0.3 


476(4) 


1507.0 


1.024(0.009) 


1.216(287) 


06 


2.03(0.02) 


0.04(0.04) 


0.292(0.004) 


0.5(0.5) 


5 +0 - 4 


471(24) 


1482.8 


1.01(0.05) 


1.242(198) 



(*) gti with biased exposure, the affected data are not taken into account in the analyses. 

The second third and forth columns report the best fit estimate for the model in Eq. [2] The fifth column reports the value of the SED peak 
analytically estimated from Eq.[2]according to the best fit results. The sixth and seventh columns report the E p and S p best fit estimates using 
as best fit model Eq.[3] In the eighth column we report the flux in the 0.3-10.0 keV band, evaluated by Xspec integrating the Eq.[2]model. The 
ninth colum reports the SED peak flux luminosity evaluated as L p = S p 4n D 2 , where D L ^ 134.1 Mpc is the luminosity distance. In the last 
column we report the reduced x 1 an d the degrees of freedom concerning the Eq.[2]fit. 
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Table 2. Swift-XRT Spectral analysis of Mrk 421. continued 
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73 
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^40C4^ 


1007 1 

1U7 / . 1 


737C0 000^ 

U. / Ji^U.UU7 ) 


QddC??l\ 

\J .y-t'-tyz.Z. 1 ) 


24 


7 o i co o n 

L.U1 ^U.U 1 1 


OQCO 04 s ! 


71 SCO 007 s ! 


SCO 7^ 


9 +0 - 2 


^44C4~1 


1077 8 

1U / Z..O 


74C0 000 s ! 

U. / t^U.UU7 J 


1 0(W??(Ti 

1 A/UJ^Z,Z.U7 


Z. J 


1 osco on 


03C0 03 s ! 


734C0 007 s ! 


7C7~1 


- -C- -"1 


- -C- -'l 

V I 


1 770 


- -C- 


1 1?7n^S^ 
i , lZj j yz. y 


26 


1.85(0.01) 


-0.02(0.03) 


0.295(0.003) 




- -(- -) 


- -(- -) 


1809.2 


- -(- -) 


1.290(269) 


27 


1.83(0.01) 


0.02(0.03) 


0.285(0.003) 


> 100 


- -(- -) 


- -(- -) 


1676.1 


- -(- -) 


1.249(268) 


28 


1.85(0.01) 


0.02(0.03) 


0.276(0.003) 


> 100 


- -(- -) 


- -(- -) 


1595.7 


- -(- -) 


1.322(261) 


29 


1.8(0.01) 


0.06(0.03) 


0.285(0.002) 


> 10 


- -(- -) 


- -(- -) 


1662.8 


- -(- -) 


1.132(293) 


30 


1.83(0.01) 


0.02(0.03) 


0.29(0.002) 


> 100 


- -(- -) 


- -(- -) 


1701.8 


- -(- -) 


1.178(284) 


31 


1.83(0.02) 


0.08(0.03) 


0.401(0.004) 


10(15) 


- -(- -) 


- -(- -) 


2260.2 


- -(- -) 


1.022(237) 


32 


1.83(0.02) 


0.08(0.03) 


0.281(0.003) 


> 10 


- -(- -) 


- -(- -) 


1598.0 


- -(- -) 


0.894(231) 


33 


1.84(0.02) 


0.11(0.04) 


0.289(0.003) 


6(5) 


6 +6 

u 2 


531(19) 


1593.0 


1.14(0.04) 


0.876(220) 


34 


1.82(0.03) 


0.17(0.06) 


0.29(0.005) 


3(2) 


3 5 +3 

J -0.8 


520(16) 


1557.4 


1.12(0.03) 


1.092(101) 



Obsld 00030352012 Date 6/20/2006 MJD 53907.001774 



01 


1.92(0.03) 


0.09(0.08) 


0.311(0.007) 


3(4) 


- -(- -) 


- -(- -) 


1629.2 


- -(- -) 


0.853(64) 


02 


1.83(0.01) 


0.16(0.02) 


0.303(0.002) 


3.4(0.9) 


■x 4+0.6 

-'•^-0.4 


538(6) 


1631.2 


1.16(0.01) 


1.179(297) 


03 


1.83(0.01) 


0.16(0.03) 


0.286(0.002) 


3(1) 


t 1+0.7 
-'-0.4 


506(7) 


1535.8 


1.09(0.01) 


1.065(247) 


04 


1.87(0.01) 


0.16(0.03) 


0.267(0.002) 


2.6(0.8) 


2 6 +0 " 5 

^•"-0.3 


455(5) 


1399.4 


0.98(0.01) 


1.113(235) 


05 


1.87(0.01) 


0.18(0.03) 


0.253(0.002) 


2.3(0.6) 


7 3+°- 3 
^-0.2 


429(5) 


1313.9 


0.92(0.01) 


1.163(227) 


06 


1.83(0.02) 


-0.01(0.03) 


0.198(0.002) 




- -(- -) 


- -(- -) 


1191.2 


- -(- -) 


0.932(218) 


07 


1.92(0.02) 


0.09(0.04) 


0.225(0.003) 


3(2) 


2 - 8 -0.7 


376(6) 


1189.5 


0.81(0.01) 


0.968(187) 


08 


1.9(0.02) 


0.14(0.05) 


0.244(0.003) 


2(1) 


9 1+0.9 
z —'-0.4 


408(7) 


1269.1 


0.88(0.01) 


1.046(149) 


09 


1.91(0.02) 


0.16(0.05) 


0.255(0.004) 


2(0.8) 


9+0.6 
-0.3 


422(7) 


1305.3 


0.91(0.01) 


1.071(112) 


10 


1.88(0.01) 


0.15(0.03) 


0.275(0.002) 


2.5(0.8) 


2 5 +as 

z " J -0.3 


464(5) 


1438.3 


1(0.01) 


1.145(239) 


11 


1.85(0.01) 


0.14(0.03) 


0.294(0.003) 


3(2) 


3 4 +1 

-''^-0.6 


514(8) 


1577.8 


1.11(0.02) 


1.284(229) 


12 


1.83(0.01) 


0.2(0.03) 


0.316(0.003) 


2.6(0.6) 


2 6 +0 - 4 

^""-0.2 


548(6) 


1657.9 


1.18(0.01) 


1.255(232) 


13 


1.87(0.01) 


0.15(0.03) 


0.325(0.003) 


2.8(0.8) 


9 a +0-5 
2 -°-0.3 


556(6) 


1713.4 


1.2(0.01) 


1.082(262) 


14 


1.85(0.01) 


0.19(0.03) 


0.337(0.003) 


2.5(0.6) 


2.5«| 


578(6) 


1762.0 


1.24(0.01) 


1.167(258) 


15 


1.91(0.01) 


0.16(0.03) 


0.319(0.003) 


1.8(0.3) 




524(5) 


1620.8 


1.13(0.01) 


1.199(252) 



(*) gti with biased exposure, the affected data are not taken into account in the analyses. 

The second third and forth columns report the best fit estimate for the model in Eq. [2] The fifth column reports the value of the SED peak 
analytically estimated from Eq.[2]according to the best fit results. The sixth and seventh columns report the E p and 5,, best fit estimates using 
as best fit model Eq.[3] In the eighth column we report the flux in the 0.3-10.0 keV band, evaluated by Xspec integrating the Eq.[2]model. The 
ninth colum reports the SED peak flux luminosity evaluated as L p - S p 4n D 2 L , where D L ^ 134.1 Mpc is the luminosity distance. In the last 
column we report the reduced x 1 and the degrees of freedom concerning the Eq.[2]fit. 



A. Tramacere,e? ah: Swift observations of the very intense flaring activity of Mrk 421 during 2006. 21 
Table 2. Swift-XRT Spectral analysis of Mrk 421. continued 



Orbit 


a b 


K 


E p * 
keV 


E P 
keV 


s P 


flux 0.3-10 keV 

10 ^ erg cm~^ s -1 


h 

10 45 erg/s 


Xf/dof 


16 


1.92(0.01) 0.14(0.03) 


0.316(0.003) 


2(0.5) 


9+0.3 
-0.2 


519(5) 


1623.1 


1.12(0.01) 


1.068(246) 



Obsld 00030352013 Date 6/22/2006 MJD 53908.050644 
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1.191(245) 


02 


1.93(0.01) 


0.15(0.03) 


0.302(0.002) 


1.8(0.3) 


1.8*« 


493(4) 


1537.5 


1.061(0.009) 


1.360(256) 


03 


1.94(0.01) 


0.18(0.03) 


0.308(0.002) 


1.5(0.2) 


1 45 +ai 


498(4) 


1522.5 


1.072(0.009) 


1.088(253) 


04 


1.98(0.01) 


0.15(0.03) 


0.269(0.002) 


1.2(0.1) 


1 2 +01 


431(4) 


1323.0 


0.927(0.009) 


1.083(241) 


05 


2.01(0.01) 


0.08(0.03) 


0.252(0.002) 


0.9(0.2) 


9 +0 2 
Kj.y _ Q3 


404(4) 


1268.0 


0.869(0.009) 


1.006(213) 


06 


1.96(0.01) 


0.12(0.03) 


0.243(0.002) 


1.5(0.3) 


1 5 +a2 
l,J -o.i 


392(4) 


1231.3 


0.843(0.009) 


1.208(235) 


07 


1.93(0.01) 


0.03(0.03) 


0.243(0.002) 


> 10 


- -(- -) 


- -(- -) 


1322.9 


- -(- -) 


1.063(245) 


08 


1.82(0.01) 


0.03(0.02) 


0.281(0.002) 


> 100 


- -(- -) 


- -(- -) 


1667.5 


- -(- -) 


1.434(277) 


09 


1.79(0.01) 


0.09(0.03) 


0.263(0.002) 


> 10 


12 + > 6 


545(32) 


1516.0 


1.17(0.07) 


1.163(265) 


10 


1.82(0.01) 


0.07(0.03) 


0.241(0.002) 


> 10 


- -(- -) 


- -(- -) 


1384.9 


- -(- -) 


1.275(256) 


11 


1.88(0.01) 


0.02(0.03) 


0.241(0.002) 


> 100 


- -(- -) 


- -(- -) 


1374.1 


- -(- -) 


1.247(256) 


12 


1.88(0.01) 


0.02(0.03) 


0.229(0.002) 


> 100 


- -(- -) 


- -(- -) 


1300.0 


- -(- -) 


1.051(244) 


13 


1.86(0.01) 


0.02(0.03) 


0.25(0.002) 


> 100 


- -(- -) 


- -(- -) 


1438.9 


- -(- -) 


1.036(255) 


14 


1.88(0.01) 


0.02(0.03) 


0.245(0.002) 


> 100 


- -(- -) 


- -(- -) 


1393.2 


- -(- -) 


1.221(254) 


15 


1.88(0.01) 


0.12(0.02) 


0.231(0.001) 


3(1) 


3 l +0 7 


395(4) 


1227.7 


0.85(0.009) 


1.334(308) 



Obsld 00030352014 Date 6/23/2006 MJD 53909.394900 



01 


1.74(0.01) 


0.02(0.02) 


0.326(0.003) 


> 100 


--(- 


-) 


--(- 


-) 


2060.1 


--(- 


-) 


1.163(281) 


02 


1.67(0.01) 


0.06(0.03) 


0.372(0.003) 


> 100 


--(- 


-) 


--(- 


-) 


2435.1 


--(- 


-) 


1.240(255) 


03 


1.66(0.01) 


0.02(0.03) 


0.411(0.004) 


> 100 


--(- 


-) 


--(- 


-) 


2783.6 


--(- 


-) 


1.110(276) 


04 


1.69(0.01) 


0.04(0.03) 


0.398(0.003) 


> 100 


--(- 


-) 


--(- 


-) 


2603.6 


--(- 


-) 


1.234(265) 


05 


1.63(0.01) 


0.08(0.03) 


0.436(0.004) 


> 100 


--(- 


-) 


--(- 


-) 


2938.2 


--(- 


-) 


1.357(284) 


06 


1.67(0.01) 


0.06(0.02) 


0.434(0.003) 


> 100 


--(- 


-) 


--(- 


-) 


2845.1 


--(- 


-) 


1.079(291) 


07 


1.68(0.01) 


0.07(0.03) 


0.444(0.004) 


> 100 


--(- 


-) 


--(- 


-) 


2870.2 


--(- 


-) 


1.045(282) 


08 


1.76(0.01) 


0.07(0.03) 


0.393(0.003) 


> 10 


--(- 


-) 


--(- 


-) 


2366.8 


--(- 


-) 


1.098(255) 


09 


1.8(0.01) 


0.03(0.02) 


0.383(0.003) 


> 100 


--(- 


-) 


--(- 


-) 


2302.1 


--(- 


-) 


1.041(286) 


10 


1.79(0.01) 


0.06(0.03) 


0.384(0.003) 


> 10 


--(- 


-) 


--(- 


-) 


2259.2 


--(- 


-) 


1.025(249) 


Ibsld 00215769000 


Date 6/23/2006 MJD 53909. 671942 


















01 


1.64(0.01) 


0.05(0.02) 


0.386(0.003) 


> 100 


--(- 


-) 


--(- 


-) 


2617.0 


--(- 


-) 


1.099(310) 



Obsld 00030352015 Date 6/24/2006 MJD 53910.068887 



01 


1.83(0.01) 


0.04(0.03) 


0.447(0.004) 


> 100 


- -(- -) 


- -(- -) 


2601.0 


- -(- -) 


1.103(252) 


02 


1.85(0.01) 


0.02(0.03) 


0.409(0.003) 


> 100 


- -(- -) 


- -(- ") 


2367.8 


- -(- -) 


1.206(263) 


03 


1.85(0.01) 


0.04(0.03) 


0.406(0.003) 


> 10 


- -(- -) 


- -(- -) 


2323.1 


- -(- -) 


1.157(263) 


04 


1.83(0.01) 


0.07(0.03) 


0.418(0.003) 


> 10 


- -(- -) 


- -(- -) 


2389.7 


- -(- -) 


1.229(259) 


05 


1.84(0.01) 


0.1(0.02) 


0.439(0.003) 


6(4) 


6* 


816(22) 


2440.4 


1.76(0.05) 


1.071(271) 


06 


1.78(0.01) 


0.12(0.03) 


0.473(0.004) 


9(7) 


9% 


957(40) 


2707.7 


2.06(0.09) 


1.454(255) 


07 


1.78(0.01) 


0.1(0.03) 


0.478(0.004) 


> 10 


12+f 


1009(61) 


2775.9 


2.2(0.1) 


1.300(252) 


08 


1.78(0.01) 


0.13(0.03) 


0.485(0.004) 


7(5) 


r* 


968(35) 


2766.7 


2.08(0.07) 


1.399(251) 


09 


1.87(0.01) 


0.1(0.03) 


0.455(0.004) 


4(3) 


4% 


800(14) 


2467.9 


1.72(0.03) 


1.054(259) 


10 


1.9(0.01) 


0.08(0.03) 


0.458(0.004) 


4(3) 


3 7 +3 

' —0.9 


780(13) 


2454.4 


1.68(0.03) 


1.033(242) 


11 


1.92(0.01) 


0.09(0.03) 


0.402(0.003) 


3(2) 


2-9!o, 


670(7) 


2128.7 


1.44(0.01) 


1.206(256) 


12 


1.93(0.01) 


0.06(0.03) 


0.376(0.003) 


4(4) 


4tl 


630(11) 


2009.5 


1.36(0.02) 


1.137(249) 


13 


1.92(0.01) 


0.07(0.03) 


0.359(0.003) 


4(4) 




609(13) 


1926.6 


1.31(0.03) 


1.065(247) 



(*) gti with biased exposure, the affected data are not taken into account in the analyses. 

The second third and forth columns report the best fit estimate for the model in Eq. [2] The fifth column reports the value of the SED peak 
analytically estimated from Eq.[2]according to the best fit results. The sixth and seventh columns report the E p and S p best fit estimates using 
as best fit model Eq.[3] In the eighth column we report the flux in the 0.3-10.0 keV band, evaluated by Xspec integrating the Eq.[2]model. The 
ninth colum reports the SED peak flux luminosity evaluated as L p - S p 4n D 2 L , where D L ^ 134.1 Mpc is the luminosity distance. In the last 
column we report the reduced x 1 and the degrees of freedom concerning the Eq.[2]fit. 



22 A. Tramacere,er ah: Swift observations of the very intense flaring activity of Mrk 421 during 2006. 

Table 2. Swift-XRT Spectral analysis of Mrk 421. continued 



Orbit a 


b 


K 


E p * 
keV 


E P 
keV 


s P 


flux 0.3-10 keV 


L„ 

10 45 ergls 




14 1.96(0.01) 

15 1.98(0.01) 

16 1.95(0.02) 


0.06(0.03) 

0.1(0.03) 

0.08(0.04) 


0.329(0.003) 
0.309(0.003) 
0.294(0.004) 


2(1) 

1.3(0.3) 

2(1) 


2 1 +1 

1 3 +02 

2 1 +2 


534(6) 
496(6) 
479(7) 


1722.2 
1567.0 
1531.0 


1.15(0.01) 
1.07(0.01) 
1.03(0.01) 


1.277(236) 
1.135(212) 
1.305(162) 


Obsld 00030352016 


Date 6/27/2006 MJD 53913.139963 












01 1.99(0.01) 

02 1.96(0.01) 


0(0.02) 
-0(0.03) 


0.253(0.002) 
0.243(0.002) 




- -(- -) 

- -(- -) 


- -(- -) 

- -(- -) 


1350.6 
1324.9 


- -(- -) 

- -(- -) 


1.001(261) 
1.079(247) 



Obsld 00219237000 Date 7/15/2006 MJD 53931.277331 



01 


1.65(0.02) 


0.14(0.04) 


0.509(0.007) 


> 10 


- -(- -) 


- <- ") 


3193.8 


- -(- -) 


1.048(171) 


02 


1.63(0.02) 


0.17(0.04) 


0.543(0.006) 


13(11) 


13 + i° 


1393(105) 


3406.2 


3(0.2) 


1.028(208) 


03 


1.66(0.02) 


0.17(0.04) 


0.534(0.006) 


9(8) 


911 


1249(80) 


3240.9 


2.7(0.2) 


0.979(191) 


04 


1.72(0.02) 


0.13(0.04) 


0.507(0.006) 


> 10 


12 +i6 


1142(90) 


3016.7 


2.5(0.2) 


1.009(191) 


05 


1.7(0.02) 


0.18(0.04) 


0.498(0.006) 


6(4) 


6 +4 

u 2 


1050(44) 


2906.0 


2.26(0.09) 


1.017(174) 


06 


1.66(0.02) 


0.13(0.03) 


0.48(0.005) 


> 10 


,0+28 


1275(134) 


3008.1 


2.7(0.3) 


0.960(228) 



(*) gti with biased exposure, the affected data are not taken into account in the analyses. 

The second third and forth columns report the best fit estimate for the model in Eq. [2] The fifth column reports the value of the SED peak 
analytically estimated from Eq. [2] according to the best fit results. The sixth and seventh columns report the E p and S best fit estimates using 
as best fit model Eq.[3] In the eighth column we report the flux in the 0.3-10.0 keV band, evaluated by Xspec integrating the Eq.[2]model. The 
ninth colum reports the SED peak flux luminosity evaluated as L p = S p 4n D\, where D L ^ 134.1 Mpc is the luminosity distance. In the last 
column we report the reduced x 1 and the degrees of freedom concerning the Eq.[2]fit. 



A. Tramacere,ef ah: Swift observations of the very intense flaring activity of Mrk 421 during 2006. 
Table 3. Swift-XRT Orbit merged spectral analysis of Mrk 421. 



23 



Interval 


a 


b 


K 


E p * 
keV 


keV 


s P 


flux 0.3-10 keV 


10 45 erg/s 


Xil&of 


Obsld 00206476000 Date 4/22/2006 


MJD 53847.252792 












01 


1.68(0.009) 


0.11(0.02) 


0.349(0.002) 


27(21) 




941(60) 


2181.3 


2(0.1) 


1.275(398) 


01 X+B 


1.68(0.01) 


0.12(0.02) 


0.351(0.002) 




20(+>°) 








1.299(401) 


02 
03 
04 
05 


1.76(0.01) 
1.8(0.01) 
1.79(0.01) 
1.804(0.009) 


0.18(0.02) 
0.09(0.02) 
0.14(0.03) 
0.13(0.02) 


0.289(0.002) 
0.26(0.002) 
0.257(0.002) 
0.281(0.002) 


5(1) 
14(14) 

5(3) 
6(2) 


4 7 +l 

^' ' -0.6 

14 + ' 5 

5 +2 
-i 

5.T_] 


558(9) 
540(27) 
488(11) 
532(9) 


1620.6 
1499.0 
1434.8 
1570.4 


1.2(0.02) 
1.16(0.06) 
1.05(0.02) 
1.14(0.02) 


1.187(300) 
1.079(331) 
1.122(254) 
1.062(344) 



Obsld 00030352005 Date 4/25/2006 MJD 53850.267940 



01 1.84(0.01) 0.15(0.02) 0.273(0.002) 3(1) 3.3^ 480(5) 1466.7 1.03(0.01) 1.350(288) 

02 1.859(0.008) 0.16(0.02) 0.266(0.001) 2.8(0.5) 2.8+°j 459(3) 1411.6 0.988(0.006) 1.363(379) 

Obsld 00030352006 Date 4/26/2006 MJD 53851.204919 

01 1.824(0.009) 0.16(0.02) 0.289(0.002) 3.6(0.9) 3.6^ 519(5) 1570.7 1.12(0.01) 1.416(336) 

Obsld 00030352007 Date 4/26/2006 MJD 53851.951963 

01 1.87(0.01) 0.16(0.03) 0.265(0.002) 2.6(0.6) 2.6^ 451(5) 1387.8 0.97(0.01) 1.085(258) 

Obsld 00030352008 Date 6/14/2006 MJD 53900.016100 

01(c) 2.05(0.01) 0.15(0.03) 0.233(0.002) 0.67(0.1) 377(3) 1102.4 0.811(0.006) 1.108(333) 

Obsld 00030352009 Date 6/15/2006 MJD 53901.490493 



01 1.86(0.01) 0.03(0.02) 0.17(0.001) > 100 --(--) --(--) 970.3 --(--) 1.191(400) 

02 1.87(0.01) 0.15(0.03) 0.191(0.001) 2.7(0.8) 2.7+g-f 326(3) 1007.0 0.701(0.006) 1.149(319) 

Obsld 00030352010 Date 6/16/2006 MJD 53902.025296 



01 


1.858(0.007) 


0.1(0.01) 


0.226(0.001) 


6(2) 


6!? 


409(6) 


1244.5 


0.88(0.01) 


1.297(413) 


02 


1.77(0.008) 


0.07(0.01) 


0.27(0.001) 


> 10 


- -(- -) 


- -(- -) 


1614.6 


- -(- -) 


1.227(420) 


03 


1.713(0.008) 


0.04(0.02) 


0.3(0.001) 


> 100 


- -(- -) 


- -(- -) 


1927.8 


- -(- -) 


1.268(423) 


04 


1.766(0.006) 


0.04(0.01) 


0.274(0.001) 


> 100 


- -(- -) 


- -(- -) 


1680.5 


- -(- -) 


1.273(477) 


05 


1.823(0.007) 


-0.02(0.01) 


0.243(0.001) 




- -(- -) 


- -(- -) 


1478.1 


- -(- -) 


1.310(439) 


06 


1.748(0.008) 


0.03(0.02) 


0.294(0.001) 


> 100 


- -(- -) 


- -(- -) 


1840.3 


- -(- -) 


1.320(414) 



Obsld 00030352011 Date 6/18/2006 MJD 53904.038766 



01 


1.808(0.009) 


0.12(0.02) 


0.319(0.002) 


7(3) 


T\ 


613(13) 


1794.7 


1.32(0.03) 


1.355(378) 


02 


1.82(0.01) 


0.19(0.02) 


0.341(0.002) 


3(0.7) 


q+0.4 
-0.3 


600(6) 


1814.2 


1.29(0.01) 


1.348(315) 


03 


1.932(0.008) 


0.12(0.02) 


0.293(0.002) 


1.9(0.3) 


I O+0.2 


479(3) 


1508.2 


1.031(0.006) 


1.331(385) 


04(c) 


1.985(0.009) 


0.09(0.02) 


0.256(0.002) 


1.2(0.2) 




410(3) 


1299.8 


0.882(0.006) 


1.115(325) 


05 


1.893(0.008) 


0.07(0.02) 


0.277(0.001) 


6(4) 




485(8) 


1506.5 


1.04(0.02) 


1.209(395) 


06 


1.884(0.006) 


0.08(0.01) 


0.271(0.001) 


5(2) 


5 +2 

-0.9 


475(5) 


1475.1 


1.02(0.01) 


1.253(473) 


07 


1.866(0.008) 


0.05(0.02) 


0.222(0.001) 


> 10 


- -(" ") 


- -(- -) 


1248.9 


- -(- -) 


1.264(419) 


08(c) 


1.989(0.008) 


0.09(0.02) 


0.213(0.001) 


1.2(0.1) 


1 2 +(u 


341(2) 


1082.0 


0.734(0.004) 


1.148(392) 


09 


1.875(0.007) 


0.01(0.01) 


0.273(0.001) 


> 100 


- -(- -) 


- -(- -) 


1569.5 


- -(- -) 


1.349(466) 


10 


1.821(0.007) 


0.06(0.01) 


0.298(0.001) 


> 10 


- -(- -) 


- -(- -) 


1720.5 


- -(- -) 


1.295(467) 



Obsld 00030352012 Date 6/20/2006 MJD 53907.001774 



01 


1.842(0.007) 


0.16(0.01) 


0.289(0.001) 


3.1(0.5) 


3 l +03 

J- -0.2 


506(3) 


1545.3 


1.089(0.006) 


1.179(419) 


02 


1.871(0.008) 


0.11(0.02) 


0.228(0.001) 


4(1) 


n 0+0.9 


398(4) 


1230.8 


0.856(0.009) 


1.238(356) 


03 


1.861(0.008) 


0.17(0.02) 


0.287(0.001) 


2.6(0.4) 


2 6 +0 - 2 

^""-0.2 


491(3) 


1509.5 


1.056(0.006) 


1.322(382) 


04 


1.857(0.009) 


0.17(0.02) 


0.331(0.002) 


2.6(0.4) 


2 6 +0 - 2 

^•"-0.2 


568(4) 


1741.4 


1.222(0.009) 


1.253(344) 


05 


1.915(0.009) 


0.15(0.02) 


0.318(0.002) 


1.9(0.3) 


I 9+0.1 


522(4) 


1623.3 


1.123(0.009) 


1.283(328) 



(c) Orbit with a strong cooling contamination. 

The second third and forth column report the best fit estimate for the model in Eq. [2] The fifth column reports the value of the SED peak 
analytically estimated from Eq.[2]according to the best fit results. The sixth and seventh columns report the E p and S ,, best fit estimates using 
as best fit model Eq.f3] In the eighth column we report the flux in the 0.3-10.0 keV band, evaluated by X-spec integrating the Eq.[2]model. The 
ninth colum reports the SED peak flux luminoisty evaluated as L p - S p 4n D 2 L , where D L ^ 134.1 Mpc is the luminosity distance. In the last 
column we report the reduced x 1 and the degrees of freedom concerning the Eq.[2]fit. 



24 A. Tramacere,er ah: Swift observations of the very intense flaring activity of Mrk 421 during 2006. 

Table 2. Swift-XRT Orbit merged spectral analysis of Mrk 421 . continued 



Interval a 


b 


K 


E p * 


E P 


s p 


flux 0.3-10 keV L p 


tf/dof 








keV 


keV 


VQT^erg cm~^ s 


Vr n ergcm~ 1 i _1 10 4s erg/s 





Obsld 00030352013 Date 6/22/2006 MJD 53908.050644 



01 


1.942(0.007) 


0.17(0.01) 


0.3(0.001) 


1.5(0.1) 


1 A7+0.O6 

1 - H 1 -0.05 


485(3) 


1493.6 


1.044(0.006) 


1.261(393) 


02(c) 


1.982(0.008) 


0.12(0.02) 


0.255(0.001) 


1.18(0.1) 


1 1 e+0.07 


408(2) 


1272.0 


0.878(0.004) 


1.247(354) 


03 


1.93(0.01) 


0.04(0.03) 


0.243(0.002) 


> 10 


- -(- -) 


- -(- -) 


1323.0 


- -(- -) 


1.119(244) 


04 


1.81(0.01) 


0.03(0.02) 


0.281(0.002) 


> 100 


- -(- -) 


- -(- -) 


1665.7 


- -(- -) 


1.403(278) 


05 


1.805(0.009) 


0.09(0.02) 


0.251(0.001) 


14(12) 


14+^ 3 


520(23) 


1446.2 


1.12(0.05) 


1.268(351) 


06 


1.87(0.006) 


0.02(0.01) 


0.242(0.001) 


> 100 


- -(- -) 


- -(- -) 


1380.5 


- -(- -) 


1.263(455) 



Obsld 00030352014 Date 6/23/2006 MJD 53909.394900 



01 


1.74(0.01) 


0.02(0.02) 


0.326(0.003) 


> 100 


- -(- -) 


--(- 


-) 


2060.1 


--(- 


-) 


1.163(281) 


02 


1.67(0.01) 


0.06(0.03) 


0.372(0.003) 


> 100 


- -(- -) 


--(- 


-) 


2435.1 


--(- 


-) 


1.240(255) 


03 


1.66(0.01) 


0.02(0.03) 


0.411(0.004) 


> 100 


- -(- -) 


--(- 


-) 


2783.6 


--(- 


-) 


1.110(276) 


04 


1.69(0.01) 


0.04(0.03) 


0.398(0.003) 


> 100 


- -(- -) 


--(- 


-) 


2603.6 


--(- 


-) 


1.234(265) 


05 


1.63(0.01) 


0.08(0.03) 


0.436(0.004) 


> 100 


- -(- -) 


--(- 


-) 


2938.2 


--(- 


-) 


1.357(284) 


05 X+B 


1.61(0.01) 


0.13(0.02) 


0.434(0.004) 




340 












1.399(287) 


06 


1.67(0.01) 


0.06(0.02) 


0.434(0.003) 


> 100 


- -(- -) 


--(- 


-) 


2845.1 


--(- 


-) 


1.079(291) 


07 


1.68(0.01) 


0.07(0.03) 


0.444(0.004) 


> 100 


- -(- -) 


--(- 


-) 


2870.2 


--(- 


-) 


1.045(282) 


08 


1.76(0.01) 


0.07(0.03) 


0.393(0.003) 


> 10 


- -(- -) 


--(- 


-) 


2366.8 


--(- 


-) 


1.098(255) 


09 


1.8(0.01) 


0.03(0.02) 


0.383(0.003) 


> 100 


- -(- -) 


--(- 


-) 


2302.1 


--(- 


-) 


1.041(286) 


10 


1.79(0.01) 


0.06(0.03) 


0.384(0.003) 


> 10 


- -(- -) 


--(- 


-) 


2259.2 


--(- 


-) 


1.025(249) 


Obsld 00215769000 Date 6/23/2006 


MJD 53909.671942 
















01 


1.64(0.01) 


0.05(0.02) 


0.386(0.003) 


> 100 


- -(- -) 


--(- 


-) 


2617.0 


--(- 


-) 


1.099(310) 



Obsld 00030352015 Date 6/24/2006 MJD 53910.068887 



01 


1.844(0.007) 


0.04(0.01) 


0.422(0.002) 


> 100 


- -(- -) 


- -(- -) 


2434.7 


- -(- -) 


1.317(418) 


02 


1.778(0.008) 


0.12(0.02) 


0.479(0.002) 


8(4) 




970(22) 


2752.4 


2.09(0.05) 


1.474(399) 


03 


1.832(0.009) 


0.09(0.02) 


0.43(0.002) 


9(7) 


Q + 6 

y -3 


829(24) 


2422.2 


1.78(0.05) 


1.181(354) 


04 


1.907(0.005) 


0.08(0.01) 


0.398(0.001) 


4(1) 


3 6 +07 

-'■"-0.4 


675(4) 


2131.8 


1.452(0.009) 


1.259(482) 


05(c) 


1.964(0.009) 


0.08(0.02) 


0.313(0.002) 


1.7(0.4) 


1 7 +02 


505(3) 


1618.7 


1.087(0.006) 


1.367(322) 



Obsld 00030352016 Date 6/27/2006 


MJD 53913.139963 












01 1.973(0.008) 0(0.02) 


0.248(0.001) - - 


- -(- -) 


- -(- -) 


1338.5 


- -(- -) 


1.072(339) 


Obsld 00219237000 Date 7/15/2006 


MJD 53931.277331 












01 1.65(0.01) 0.17(0.02) 
01 X+B 1.64(0.01) 0.20(0.02) 


0.532(0.004) 11(5) 
0.535(0.05) 


U + i 
8C?) 


1299(50) 


3288.2 


2.8(0.1) 


1.147(322) 
1.235(325) 


02 1.684(0.009) 0.16(0.02) 


0.504(0.003) 9(4) 


9tl 


1150(32) 


3035.6 


2.47(0.07) 


1.136(369) 



(c) Orbit with a strong cooling contamination. 

The second third and forth column report the best fit estimate for the model in Eq. [2] The fifth column reports the value of the SED peak 
analytically estimated from Eq. [2] according to the best fit results. The sixth and seventh columns report the E p and S p best fit estimates using 
as best fit model Eq.[3] In the eighth column we report the flux in the 0.3-10.0 keV band, evaluated by X-spec integrating the Eq.[2]model. The 
ninth colum reports the SED peak flux luminoisty evaluated as L p = 5 ; , 4n D 2 L , where D L ^ 134.1 Mpc is the luminosity distance. In the last 
column we report the reduced x 1 and the degrees of freedom concerning the Eq.[2]fit. 
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